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ABSTRACT 


NOAA-6 satellite AVHRR data and AXBT data were collected 
in the Northeast Pacific Ocean in late 1980 as part of the 
Naval Postgraduate School-sponsored Acoustic Storm Transfer 
and Response Experiment which was in turn part of the U.S.- 
Canadian Storm Transfer and Response Experiment (STREX). Some 
of the problems in transferring AXBT geographical positions to 
satellite images were solved by designing a computer program 
with accuracies of less than 2 pixels. Thermal comparisons 
were made between AXBT, NOAA-6, and GOSSTCOMP data with the 
result that NOAA-6 data was on the average 2.9?C colder than 
AXBT data and 3.2?C colder than GOSSTCOMP data. Linear regression 
methods reduced to 0.39С the difference between NOAA-6 and AXBT 
data. Use of this method over a period of 15 days produced a 
En error of 0.5°C. 

Although NOAA-6 cannot sense directly the subsurface thermal 
Structure, it is excellent for observing surface manifestations 
of horizontal thermal features. Further investigation into 
using satellite data as the basis of an empirical relationship 
between the surface temperature and the subsurface vertical 


thermal structure is warranted. 





TABLE OF CONTENTS 


INTRODUCTION -------------------------------------‏ ا 


D. КОО ERENS TNC TIN OCEANOGRAPHIC RESEARCH --------- 


A. 


OCEANOGRAPHIC CHARACTERISTICS OF THE 
PROJECT AREA --------------------------------- 


l.  Thermal-Salinity Structure --------------- 
DI SNS ————-——-—-—-——-——--------- 
3 Eronts and Eddies ------------------------ 
ШЕЕК ОРИ een Teo IN OCEAN THERMAL STUDIES === 


1. Problems Associated with a Satellite 
Ра+а Вазе -------------------------------- 


a. Atmospheric Attenuation -------------- 
b. Location Accuracy -------------------- 


C. Remote Sensing of the Vertical 
Structure ---------------------------- 


NOAA-6 OPERATION ----------------------------- 
l. The Spacecraft --------------------------- 
a. Physical Structure ------------------- 


b. The Attitude Determination and 
Control Subsystem (ADACS) ------------ 


C. Data Handling Subsystem ————— لد ت‎ 
ТТТ ІС Sens os =================== 


3. | NOAA-6 Orbital Parameters ---------------- 


III. DATA COLLECTION AND PROCESSING TECHNIQUES -------- 


A. 


AXBT COLLECTION AND PROCESSING --------------- 


l. The Air-Dropped Expendable 
Bathythermograph ------------------------- 


2. AXBT Data Processing --------------------- 


Ty 


19 


20 


22 


23 


30 





PV. 


В. SATELLITE DATA SET SELECTION AND PROCESSING --- 22 
l. NOAA-NESS --------------------------------- 55 
он o o e a 58 
b. Ephemeris Data Set -------------------- 59 
C. Initial Satellite Pass Selection ------ 64 
2.  NASA-Ames Research Center ----------------- 65 
a. IDIMS --------------------------------- 66 
b. Landmark Identification --------------- 67 
(1) Count Values --------------------- 68 
C. Pixel Identification ------------------ 69 
3. Satellite Image Navigation ---------------- "as 
а. Zoom Transfer Scope ------------------- 17 
b. IDIM's TRNSFORM ----------------------- 78 
С. NOAA-6 Orbital Dynamics --------------- 79 
d. Computer Navigation Program ----------- 83 
(1) Preliminary Programs ------------- a4 
2) Common Case Geometry ------------- 84 
(3) Units and Notation --------------- 87 
(4) Nodal and Subsatellite Point 
Calculations --------------------- 88 
(5) Buoy Pixel Identification -------- 09 
С. GOSSTCOMP ------------------------------------- ИІ 
RESULTS ------------------------------------------- 77 
A. NAVIGATION ACCURACY --------------------------- з 
B. THERMAL COMPARISONS --------------------------- 3(9 
l. Horizontal Distribution ------------------- KEB 
2. Vertical Distribution --------------------- 01 


CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY - 142 








APPENDIX A: 


Count-to-Temperature Conversion Table ----- 147 


APPENDIX B:  SCANLINE Computer Program ----------------- 149 

APPENDIX C:  TAPEDUMP Computer Program ----------------- LoL 

APPENDIX D:  AREAMAP Computer Program ------------------ 156 

APPENDIX E: LOCATE Computer Program -- ----------------- TEL 

EMEN s se M == ———————— = ۲81 

ОИ E ЕТЕР -----------........-......-.......- 187 
1 





Table 


10 


iol OF TABLES 


NOAA-6 AVHRR Channelization --------------------- 44 
NOAA-6 AVHRR Channel 4 Nonlinearity Errors ------ 47 
Біеісй-5сасіоп Ғоспас ---------------------------- 58 
Selected Satellite Passes ----------------------- 66 
БИ ЕТТЕР ۲۹۲۸٥66  × 2 > ا سے سے سے سس اہ تب یہ‎ Sul 
Satellite Data Set Orbital Parameters ----------- 82 
NAS ISO SIS ----------...................-. 87 
Statistical Summary of Navigation Accuracy ------ I 


Navigation Errors Associated with a 
ي ي ي ي ي‎ o ae nn nn 114 


Temperature Comparison Statistics --------------- 1219, 





FIGURE 


10 
ЦЫ 
12%; 
АЗ 
14 


> 


16 
7 
ШЗ 
19 
20 


21 


LISE OE FIGURES 


Major surface currents and water masses of 
the Northeast Pacific Ocean --------------------- 


Seasonal oceanographic structure of the 
Northeast Pacific Ocean ------------------------- 


November mean vertical temperature structure 
in the project area ----------------------------- 


December mean vertical temperature structure 
in the project area ----------------------------- 


Annual surface salinity of the project area ----- 
November mean layer depth in project area ------- 
December mean layer depth in project area ------- 
Atlas-F launch vehicle -------------------------- 
NOAA-6 spacecraft ------------------------------- 
NOAA-6 data flow diagram ------------------------ 
NOAA-6 AVHRR channel 4 spectral response curve -- 
AXBT patterns for the project area -------------- 
Digitizer AXBT profiles, an example ------------- 
Digitizer depth-of-isotherm summary, an example - 


Digitizer temperature-at-5-meter-intervals 
summary, an example ----------------------------- 


Header record--field-station format ------------- 
Data record--field-station format --------------- 
Ephemeris data set, an example ------------------ 
Ephemeris data set, an example ------------------ 
NOAA-6 ascending pass, from IDIMS processing ---- 


NOAA-6 descending pass, from IDIMS processing --- 


2l 


24 


25 


26 


2 


28 


29 


38 


40 


42 


46 


49 


23 


54 


56 
60 
61 
62 
63 
72 


Dis 





FIGURE 


22 


23 


24 
25 
26 
27 
28 
29 
20 
81 
12 
8 3 
34 
55 


36 


1ذ 


38 


B9 


40 
41 
42 
43 


44 


LCP TE ОЦПЕОШС ОЕ сошп values, San Francisco 
Bay entrance; cloud-free image ----------------- 


PICERINT output of count values, San Francisco 
Bay entrance; cloud-covered image -------------- 


NOAA-6 orbital plane inclination --------------- 
Satellite altitude determination --------------- 
Determination of great circle distances -------- 
Case l--orbital characteristics ---------------- 
Case 2--orbital characteristics ---------------- 
Case 3--orbital characteristics ---------------- 
Case 4--orbital characteristics ---------------- 
Case l--pixel determination -------------------- 
Case 2--pixel determination -------------------- 
Case 3--pixel determination -------------------- 
Case 4--pixel determination -------------------- 
Box geometry ----------------------------------- 


Sea surface temperature comparisons, 17 
November 1980, center track -------------------- 


Sea surface temperature comparisons, 1 
December 1980, center track -------------------- 


Sea surface temperature comparisons, 5 
December 1980, center track -------------------- 


Sea surface temperature comparisons, 95 
December 1980, northern track ------------------ 


Meander, 17 November 1980 ---------------------- 
Surface weather chart, 17 November 1980 -------- 
Surface weather chart, l December 1980 --------- 
Surface weather chart, 5 December 1980 --------- 
Vertical thermal structure, 17 November 


1980, center track ----------------------------- 


l0 





FIGURE 


45 


46 


47 


Vertical thermal structure, l1 December 1980, 
Se ae ae i c 0 —_-_---------------- 137 


Vertical thermal structure, 5 December 1980, 
د‎ = 1039 


Vertical Thermal structure, 5 December 1980, 
En uM ee ee ee سد سد‎ S 


LL 





ADACS 


ASCII 


ASTREX 
AVHRR 


AXBT 


CDA 
cm 
DCS 
DTG 
DIR 
ESA 
ESL 
GMT 
GOES 


GOSSTCOMP 


HEPAD 
ВЕЕ 2 
ШЕРТ 
DIMS 
ШОУ 
EMP 


IR 


TABLE OF ABBREVIATIONS 


Attitude Determination and Control Subsystem 


American National Standard Code for 
Information Interchange 


Acoustic Storm Transfer and Response Experiment 
Advanced Very High Resolution Radiometer 
Air-Dropped Expendable Bathythermograph 

Celsius 

Command and Data Acquisition 

centimeters 

Data Collection System 

Date Time Group; day of month followed by GMT 
Digital Tape Recorder 

Earth Sensor Assembly 

Electromagnetic Systems Laboratory 

Greenwich Mean Time (also known as Zulu time) 
Geostationary Operational Environmental Satellite 


Global Operational Sea Surface Temperature 
Computation 


High Energy Proton and Alpha Detector 

High Resolution Infrared Radiation Sounder 
High Resolution Picture Transmission 
Interactive Digital Image Manipulation System 
Instantaneous Field-of-View 

Instrument Mounting Platofrm 


Infrared 


2 





A 


MEPED 


MIERP 


MSU 


NASA 


NEAT 


NESS 


NL 


nm 


NOAA 


NS 


RSS 


SEM 


SSU 


STREX 


Шер 


WIP 


TOVS 


VHRR 


XSU 


ып 


Kelvin 

kilometer 

Meter 

Medium Energy Proton and Electron Detector 
Manipulated Information Rate Processor 
Microwave Sounding Unit 

National Aeronautics and Space Administration 
Noise Equivalent Differential Temperature 
National Environmental Satellite Service 
Scan line number 

nautical mile 

National Oceanic and Atmospheric Administration 
Sample number 

Reaction Support Structure 

sample standard deviation 

Space Environment Monitor 

Stratospheric Sounding Unit 

Storm Transfer and Response Experiment 
Total Energy Detector 

TIROS Information Processor 

TIROS Operational Vertical Sounder 

Very High Resolution Radiometer 
Cross-strap Unit 


microns 


152 





TABLE OF SYMBOLS 


ENGLISH SYMBOLS 


a semi-major axis 

an ascending node 

р semi-minor axis 

e eccentricity 

E non-rotating (fixed) earth measurement 
H mean satellite altitude 

ï inclination 

- retrograde inclination 
B1,2,3,4 box latitudes 

E AXBT (buoy) latitude 

а check procedure latitudes 
Lo Landmark latitude 

E pixel latitude 

Le subsatellite point latitude 
NL scan line number 

NS sample number 

R local earth radius 

t time 


EEK SYMBOLS 
а common angle 
Y common angle 


change in longitude between the landmark and the 
subsatellite point 


Ez 
>” 


change in the ascending node longitude‏ نت 


14 





En 


AA 
р 


E. 


change in the descending node longitude 


change in longitude between pixel and subsatellite 


AIDE 


change in longitude due to earth's rotation 


common angle 

geocentric angle 

zenith angle 

scan angle 

box longitudes 

ascending node longitude 

fixed-earth ascending node longitude 
AXBT longitude 

descending node longitude 

fixed-earth descending node longitude 
landmark longitude 

pixel longitude 

subsatellite point longitude 

great circle distance 

great circle distance between node and 


great circle distance between node and 
subsatellite point 


great circle distance between node and 
subsatellite point 


IS 


landmark 


the 





ACKNOWLEDGEMENTS 


The opportunity to work with experts in fields that had been 
only an avocation prior to my arrival at the Naval Postgradu- 
ate School has been a rewarding and enjoyable experience. 

The guidance and support of my thesis advisor, Dr. Glenn Jung, 
Professor, Department of Oceanography, and LCDR Cal Dunlap, 
Assistant Professor, Department of Oceanography, has been 
tremendous, especially when the proverbial "light at the end 
of the tunnel" was pretty dim. 

This thesis presented many problems that may have been left 
unresolved were it not for the expertise and understanding of 
Dr. James Mueller, Adjunct Professor, Department of Oceanography. 
His conceptualizations of the overall geometric problems in 
correcting satellite imagery were extremely nelpful. In addi- 
tion, satellite support came from Mr. Bob Wriglev, Technology 
Applications Branch at NASA Ames Research Center, who gave un- 
selfishly of his time both to ensure access to the NASA-1DIMS 
system and to coach me in its use. Mr. Larry Breaker, NOAA- 
NESS, was kind enough to supply the satellite tapes and the 
ephemeris data in addition to answering my many questions on 
NOAA-6 operations. 

Finally, a special thanks goes to my colleagues in the 
Environmental Acoustic Research Group, LCDR Thom Holt and 
Mr. Gene Brown (NAVOCZANO), who not only arranged the P-3C 
۰۰۰۰٦٢٣۶۹ but also provided a sounding board on which to bounce 


ldeas. 


156 





E INIRODUCTION 


The collection of oceanographic data has always been time 
consuming and expensive. With the development of environmental 
satellites, a method for the rapid gathering of oceanographic 
data was available to complement that data gathered on re- 
search cruises. Unfortunately, there are still some problems 
in using satellite data, such as the effects of the atmosphere 
on the radiative transmission path and the effect of the geo- 
metric distortion of thermal features found on polar-orbiting 
satellite images as from NOAA-6. 

This thesis is part of a series of on-going studies at the 
Naval Postgraduate School by the Department of Oceanography 
Environmental Acoustic Research Group. The overall goal of 
the Group is to continue the development of those aspects of 
acoustical oceanography that have a significant effect on naval 
tactical applications. In pursuit of this goal, the Group was 
۲٣ rticipant in the joint U.S.-Canadian Storm Transfer and 
Response Experiment (STREX) held in the fall of 1980 in the 
Memeneast Pacific Ocean. The Group is particularly interested 
in investigating whether or not satellites can fulfill the 
role presently played by ships and/or air-dropped expendable 
bathythermographs (AXBT) in gathering sea surface temperature 
data for use in forecasting the ocean's thermal structure, 
using procedures similar to those developed by the U.S. Naval 


Oceanographic Office Antisubmarine Warfare Environmental 


ШУ 





Prediction (ASWEPS) Program in the early 1960's. If direct 
correlations could be found between satellite-derived sea sur- 
face temperatures and the vertical thermal structure, then a 
rapid method of surveying the world's oceans could result, 
with numerous naval ramifications. 

As a part of the experiment conducted by the Group under 
the title Acoustic Storm Transfer and Response Experiment 
(ASTREX), this thesis was directed toward the examination of 
some of the problems in using satellites to observe the hori- 
zontal and vertical thermal distributions in the waters of 
the Northeast Pacific Ocean. Of particular interest was the 
problem of locating open-ocean geographic positions (ship, 
AXBT, etc.) on satellite images with as much accuracy as 
possible. A major portion of this thesis is devoted to this 
subject. If comparisons are to be made between satellite- 
obtained thermal values and ground-truth thermal values, then 
an elimination of location errors between the two media makes 
the results that much more significant. The reasons why NOAA- 
6 satellite imagery is distorted and a method to eliminate any 
location errors successfully are presented below. In addi- 
tion, once location accuracy was assured, various meaningful 
comparisons were made between satellite, bathythermograph, 
ENSESGOSSTCOMP (Global Operational Sea Surface Temperature 
Computation) data. The results of these comparisons are 


also presented. 
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II. REMOTE SENSING IN OCEANOGRAPHIC RESEARCH 


In 1870 and 1879 respectively, the authors Edward Everett 
Hale and Jules Verne wrote about placing an artificial satellite 
in orbit, Hale using a huge waterpowered flywheel and Verne 
a gun of sufficient muzzle velocity (Corliss, 1967). Hale 
envisioned the satellite as an aid to both navigation and 
 апасаслоп. Up until 1935, the idea of launching man into 
space remained the ideas of small amateur organizations with 
some notable exceptions, such as the efforts of Robert Goddard 
in Auburn, Massachusetts, in the 1920's. 

With the stirrings of war in Europe in the late 1930's, 
the now infamous V-2 rockets were developed and launched from 
Peenemuende by a team of German scientists including Wernher 
von Braun.  Captured en masse in 1945 by the Allies, this group 
of scientists and their hardware were transferred to the United 
States. The Army Air Force subsequently commissioned a study 
by Rand Corporation, who reported in 1946 for the first time 
Ehat an earth-orbiting artificial satellite could be used 
scientifically in the fields of meteorology, biology, and com- 
munications (Corliss, 1967). This launched a nine-year effort 
by lobbyists both inside and outside the government culminating 
in President Eisenhower's announcement on 29 July 1955 that 
the U.S. would launch an earth-orbiting scientific satellite 
to investigate the environment. The Soviet Union announced 


Similar plans on the following day. October 4, 1957 saw the 
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ШЕШЕН от Sputnik-I by Russia followed on 31 January 1958 by 
the launch of the first U.S. satellite, Explorer-I. Since 
that time, hundreds of scientific satellites have been launched 
in order to investigate topics from numerous fields. One of 
these, NOAA-6, was the satellite whose data were used on this 
project.  NOAA-6 is an earth-orbiting, environment-sensing 
spacecraft with applications in meteorology and oceanography. 
This section reviews the use of satellites in oceanographic 
studies, provides a description of the operation of NOAA-6, 
and briefly summarizes the normal oceanographic conditions for 
that region of the Northeast Pacific Ocean where the experi- 


ment took place. 


A. OCEANOGRAPHIC CHARACTERISTICS OF THE PROJECT AREA 

The area chosen for this project encompassed that region 
of the Northeast Pacific Ocean between latitudes 40 N and 50 N 
and between longitudes 126 W and 139 И. See Figure l. The 
oceanographic conditions of this region have been extensively 
ШЫ ед by Tully (1961; 1964), Tabata (1964; 1965; 1978), and 
Roden (1975) among others. The reader is referred to these 
works for more detailed information as only a brief summary 
of their findings is described below. 

As seen in Figure l, the project area is located mainly 
in an oceanic water mass transitional region between the Sub- 
arctic Water Mass, predominantly to the north of 45 N, and 
the Pacific Equatorial Water Mass, predominantly to the south 


ОЕ 23 М. The Subarctic Current flowing eastward along 45 N 
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Figure l. Major surface currents and water masses of 
the Northeast Pacific Ocean (after Kibblewhite 
3660.1900 
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latitude divides on the northwest side of the project area 
into the Alaska Current, which circles counterclockwise to 
the north along the Canadian and Alaska coastlines; and into 
the California Current, which flows southward along the western 
coast of the United States. Observations from the project area 
would be expected to show the physical characteristics of the 
Subarctic Water Mass; however, those observations in the 
southern portion of the project area may be somewhat tempered 
by the colder offshore waters of the California Current. 

ШІ Thermal-Salinity Structure 

An excess of precipitation over evaporation of approxi- 
mately 25 cm/year (Tabata and Giovando, 1962) has helped to 
create a layer of water extending to a depth of 100 meters in 
the Subarctic Water Mass that is isohaline during the winter 
months. A permanent halocline extending from 100 to between 
200 and 300 meters exists in which the salinity increases by 
EN to approximately 33.8 ?/4, and marks the maximum limit 
of seasonal effects (Tully, 1964). 

The top of the permanent halocline at 100 meters also 
marks the maximum depth of the seasonal thermocline. During 
the summer, the thermocline forms between 25 and 50 meters, 
influenced heavily by wind mixing effects alone. With the 
coming of the fall and winter months with their intensive 
storms, the surface waters begin to cool considerably and both 
convecticn and wind mixing erode the thermocline until iso- 
thermal conditions exist to the top of the permanent halocline. 


This condition usually is reached by February at which point 
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the waters continue to cool until the end of March, when the 
heating season begins. See Figure 2 for a general depiction 
of the seasonal structure. Figures 3 and 4 show the expected 
mean thermal structure for the project area for the months of 
November and December. Figure 5 shows the annual surface 
salinity maxima while Figures 6 and 7 show the expected layer 
depths also for November and December. 

Both the Subarctic Current and the California Current 
have surface speeds less than one knot with volume transport 
averaging 10 to 15 million cubic meters per second (Knauss, 
1978). As a result, the water in the project area is exposed 
to constant climatic conditions over many months and has suffi- 
cient time to adjust to seasonal variations. During November 
and December, there is an expected net heat loss in the pro- 
ject area ranging up to 400 g-cal/cm^/day (Tabata, 1961). This 
heat loss is directly responsible for the winter convective 
mixing process mentioned above; therefore, the vertical thermal 
structure is due more to the area's heat budget, storm cycle, 
and salinity layers than to any influx of new water. 

2. Internal Waves 

Internal waves have been shown to cause amplitude 
fluctuations of 4.5 to 5 meters at the level of the seasonal 
thermocline in the Northeast Pacific Ocean (Tabata and Giovando, 
1962). These oscillations may cause a periodic thickening of 
the thermocline from 10 to 50 meters due to phase differentials 
between the top and bottom of the thermocline (Tully, 1964). 


These periodic fluctuations vary with depth as oscillations of 
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Figure 2. 
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Seasonal oceanographic structure of the 
Northeast Pacific Ocean (from Tully, 
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Figure 6. November mean layer depth in project 
area (from Robinson, 1976) 


depths in feet 





28 





ure 1< December mean layer depth in project 
area (from Robinson, 1976) 


depths in feet 





23 





isotherms at the level of the halocline are five times 
greater than oscillations in the region of the thermocline. 
3. Fronts and Eddies 

The Subarctic Front, usually found between latitudes 
40 N and 45 N, may be present in the center of this project 
area. This front is characterized by the lack of a density 
front in the upper 100 meters, by the region of the strongest 
surface baroclinic flow being to the south of the surface 
temperature and salinity fronts, and by the mixed layer depth 
extending to the top of the halocline at 100 meters on its 
northern edge (Roden, 1975). In the southern section of the 
project area, eddy formation may be present in that area in- 
fluenced by the California Current which has been described 
as a relatively shallow meandering current with alternating 


EE and cold tongues (Bernstein et al., 1977). 


ВЕ OF SATELLITES IN OCEAN THERMAL STUDIES 

In 1968, a study was done comparing satellite-obtained sea 
surface temperatures with monthly mean surface temperatures 
with the result that the satellite values were anywhere from 
3 to 8.3 degrees C lower than the mean (LaViolette and Chabot, 
1968). The relative horizontal gradients observed in their 
satellite data, however, were fairly consistent with similar 
mean gradients from the historical data. This pattern of 
satellite-obtained sea surface temperatures being lower than 
the mean or the actually observed sea surface temperatures 


БЕ 516555 until today, except that technological advances have 


30 








reduced the differences in temperature between the two sets 
of data so they now are between 0.5 to 3.0 degrees C (Rao et 
шаша? Brower et al., 1975; McMillin, 1975; Cogan and 
Willard, 1976; Barnett EE Lov, Labata and Gower, 1980). 

Prior to 1972, the oceanographic use of satellite-obtained 
sea surface temperature was severely limited by both the en- 
gineering characteristics of the satellite radiometers and by 
the environmental aspects causing atmospheric attenuation. 
Large instantaneous fields-of-view (IFOV) limited the resolu- 
tion capacity of the satellite and the large values for the 
variations in the electronic signal (NEAT) caused spatial and 
temporal errors, making it difficult to detect the gradients 
associated with oceanic fronts (Legeckis, 1978). 

Several methods were proposed to remove the atmospheric 
contamination responsible for the majority of the difference 
between satellite and observed sea surface temperature values. 
The 3 to 8.3 degree difference found by LaViolette and Chabot 
(1968) came from satellite data that were not corrected for 
atmospheric attenuation, but in 1969 they developed a daily 
averaging method to lessen its impact (LaViolette and Chabot, 
mee. Vukovich (1971) developed a filtering technique to 
Seeomplish the same purpose while Smith et al., (1970) used 
a statistical method which, when compared with ship observa- 
tions, had both bias and random errors of less than 1 degree 
EEUSEng early NIMBUS satellite data. Maul and Sidran (1973) 
investigated the effects of the atmosphere, nadir angle, cloud 


amount, cloud height, and random noise which resulted in a 
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theoretical error (2 degrees C) for the NOAA satellite series, 
then soon to be launched. 

In early 1970, NOAA launched ITOS-1 which was the first 
satellite in the NOAA series of satellites of which both NOAA- 
6 and NOAA-7 are now in orbit. In early 1970, NOAA-NESS be- 
gan working on a satellite data processing model, to include 
the effects of atmospheric attenuation, which was the prede- 
cessor to the GOSSTCOMP (Global Operational Sea Surface Tempera- 
Eure Computation) model (Brower et al., 1976). With the 
launch of NOAA-2 in 1972, a more advanced radiometer was put 
into use with an IFOV of about 1 kilometer and a much reduced 
system NEAT of less than 3.0 degrees C (Legeckis, 1978). With 
this improved system, sea surface temperature fronts could be 
detected and monitored. Among the studies done during the 
following few years were those of LaViolette (1974) on upwell- 
ings off the west coast of Africa, Stumpf and Rao (1975) on 
tracking eddies in the Gulf Stream, and Bernstein et al., (1977) 
on the comparison of eddies in the California Current with 
direct observations. 

By the mid-1970's NOAA-NESS had refined their satellite 
data processing model; however, comparisons with observed data 
NOAA itself and by others found that the quality of measure- 
ments varied with time and geographical area and were related 
to the temperature gradient field; good correlation Came from 
regions of weak gradients and marginal results came from regions 
Of strong gradients (Brower et al., 1976). Klein (1979) found 


that NOAA-5 sea surface temperatures in the Northeast Pacific 
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Ocean that had been subjected to the GOSSTCOMP model were 
biased 3.5 to 3.9 degrees C and suggested that the error was 
a result of overcorrection by the model for atmospheric attenu- 
ation. With the launching of TIROS-N in 1978, NOAA-NESS up- 
dated GOSSTCOMP to take advantage of the Advanced Very High 
Resolution Radiometer (AVHRR) on this, and on the follow-on 
NOAA-6 and NOAA-7, satellites. Whereas NOAA-5 had a NEAT of 
1 to 1.5 degrees C, the TIROS-N/NOAA A-G satellite series has 
a NEAT of 0.12 degrees C (Schwalb, 1978). The improvement in 
NEAT should result in a better correlation between observed 
and satellite-derived sea surface temperatures. Chahine (1980) 
Suggested that an absolute accuracy of 1 degree C in these 
differences could be obtained by simultaneous observations of 
atmospheric and surface emissions with multi-channel radiometers, 
uSing spectral regions of the 3.7 um carbon dioxide windows as 
the main sounding channel. An instrument to accomplish this 
has yet to fly on a satellite. 
1. Problems Associated with a Satellite Data Base 

Briefly described below are three common problems asso- 
Stated with using satellites in thermal studies. Atmospheric 
attenuation is important when interpreting satellite-derived 
temperatures, while location accuracy is important when ther- 
mal comparisons are made between ship, satellite, and AXBT 
data. The depth to which present-day radiometers sense the 


thermal structure concludes the section. 
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a. Atmospheric Attenuation 

Ae described in detail in Section II.C.2, 
data from the 10.5 to 11.5 um infrared channel on NOAA-6 were 
used on this project. Radiation in this spectral region 
emitted from the earth's surface or from cloud tops is attenu- 
ated in its passage through the atmosphere to the radiometer. 
The major contribution to this attenuation is water vapor which 
can be responsible for up to a 9.0 degree C correction in the 
ЕЕ ute data (Brower et al., 1976). The amount of water 
vapor in the atmosphere varies horizontally, vertically, and 
in time with the least amount of absorption around the 9.5 to 
10.5 ym region (Fett and Mitchell, 1977). Other absorbers and 
their possible corrections are carbon dioxide (0.1 to 0.2 
degrees), ozone (0.1 degrees), and aerosols (0.1 to 0.95 degrees). 
Details on the physics of this absorption process can be found 
EN berts et al., (1976) and Weinreb and Hill (1980). 

Many atmospheric correction techniques have been 
tried in an attempt to correct satellite data. Some of these 
were discussed previously. A knowledge of the vertical mois- 
ture field would help significantly in reducing the attenuation 
effects but these data are not generally available. In any 
case, the multispectral approach to this problem seems to offer 
the best chance to reduce this type of error significantly 
(Chahine, 1980; Deschamps and Phulpin, 1980). 

"XE осал оп Accuracy 
Aa or portion of this project was devoted to 


locating geographic positions correctly on satellite imagery. 
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Estimates of location error vary widely. The data archived 
from early satellites in the NOAA series had accuracies within 
20 kilometers along the orbital track but for high zenith 
angles the accuracy decreased to within 40 kilometers (Conlon, 
1973). Subpixel accuracies were discussed by Bernstein and 
Ferneyhough (1975) on LANDSAT imagery. A satellite image con- 
taining land is usually much easier to correct geometrically 
than an image that contains open ocean, simply because control 
points are much easier to identify on the land image. 

A method of correcting VHRR imagery using a simple 
algorithm by Legeckis and Pritchard (1976) had a mean accuracy 
of 5 kilometers. A technique whereby ship positions were 
transferred to satellite images by comparing cloud features 
yielded errors of 80 to 90 kilometers although this was not the 
main purpose of the study (Cogan and Willand, 1976). Another 
study, using both a zoom transfer scope and NOAA-supplied photo- 
graphic prints of satellite imagery used in conjunction with 
plastic overlays containing latitude and longitude lines, had 
a location accuracy to within 10 to 50 kilometers depending on 
the distance of the feature being located from land control 
points (Tabata and Gower, 1980). 

pa discussion of location accuracy can be found 
in Section IV below. 

C. Remote Sensing of the Vertical Structure 

The sea "surface" temperature that a satellite 

senses is a manifestation of a complex process that occurs in 


the top few millimeters of water called the thermal boundary 
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layer. This layer is subject to the processes of net upward 
heat flux, infrared and solar radiation, and turbulence with 
the resulting temperature difference between the top and bottom 
of this layer of up to 1.0 degrees C (Katsaros, 1980). Typical 
radiometers sense only the radiation emitted from a depth of 
SU 50 um. 

Direct measurement by satellites of the deeper 
vertical thermal structure is not possible with the instruments 
carried onboard the satellites in orbit today. Techniques 
using Raman lidar systems have been developed theoretically 
and prototypes experimentally tested with reported accuracies 
within 0.2 degrees C (theoretical best value) to depths of 30 
meters (Leonard et al., 1979). Conclusions from this study 
suggest that the structure to depths of 100 meters may be de- 
tectable. The physics of the Raman spectra used in this pro- 


cess can be found in Murphy and Bernstein (1972). 


C. NOAA-6 OPERATION 

The NOAA-6 satellite is the second satellite in a series 
of third generation, polar-orbiting satellites that began with 
the launch of TIROS-N on 13 October 1978 at the Air Force 
Western Test Range, Vandenberg Air Force Base, California. The 
TIROS-N/NOAA A-G satellite series, of which NOAA-A became re- 
designated NOAA-6 upon its successful launch, is a joint re- 
search effort of the United States, the United Kingdom, and 
France and is operated by the National Environmental Satellite 


Service of the National Oceanic and Atmospheric Administration 
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(NOAA-NESS) under the U.S. Department of Commerce. The United 
Kingdom provided one of the three sounding units onboard the 
satellite, France supplied the onboard data collection system 
(DCS), the National Aeronautics and Space Administration (NASA) 
funded the development and launch of TIROS-N, and NOAA supplied 
the funds for the NOAA-6 satellite. The mission objective of 
this satellite series that directly relates to this thesis is 
the continuous monitoring of the environmental features in the 
western hemisphere which is accomplished in conjunction with 

a Second satellite system, also operated by NOAA, the Geo- 
stationary Operational Environmental Satellite (GOES) System. 
It should be noted that TIROS-N ceased operation in late 1980. 
NOAA-6 was still functioning at the writing of this thesis and 
NOAA-7 began operating in June 1981. 

For the purposes of this project, only those spacecraft 
Systems that were extensively used or are important to the 
understanding of the results are explained below. The reader 
BENE crred to Schwalb (1978), Hussey (1979), Lauritson, et al., 
(1979), and ITT Aerospace (undated) for a fully detailed des- 
cription of the many instruments onboard NOAA-6. Sections of 
these references, especially the works of Schwalb and Hussey, 
were used extensively below. 

l. The Spacecraft 

NOAA-6 used an Atlas-F launch vehicle which is a com- 
paratively small rocket approximately 28 meters tall and 
weighing about 600,000 kilograms. See Figure 8. The main body 


of the rocket detaches after launch and a second stage solid 
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Figure 8. Atlas-F launch vehicle (from Hussey, 1979) 
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rocket motor, an integral part of the NOAA-6 satellite itself, 
шый until depletion putting the satellite into a nominal 
833-kilometer orbit. 
a. Physical Structure 

The satellite itself, as shown in Figure 9, con- 
sists of three sections. The Reaction Support Structure (RSS) 
includes the injection motor mentioned above, the attitude 
control propulsion system, and an 11.6 square-meter solar cell 
array. The Instrument Mounting Platform (IMP) includes the 
attitude control sensors and the Advanced Very High Resolution 
Radiometer (AVHRR). The five-sided central structure, located 
between the RSS and the IMP, includes twelve thermal control 
louvres and the earth-facing communications antennae. The 
satellite is 3.71 meters long and 1.88 meters in diameter. Its 
weight at launch was 1420 kilograms which reduced to 737 kilo- 
grams once established in its orbit. 


b. The Attitude Determination and Control Subsystem 
(ADACS) 


When a satellite sensor, such as the AVHRR, scans 
the surface of the earth, the attitude of the spacecraft is 
extremely important in determining during data analysis just 
where the sensor looked. Any roll, pitch, or vaw on the sat- 
ellite will make the application of scan geometry extremely 
difficult and significant errors would result. Because of 
this, the ADACS system was designed to maintain the attitude 
of the spacecraft to within 0.2 degrees (3-sigma) of the local 


geographic reference (Schwalb, 1978). This value is obtained 
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through the use of three mutually-orthogonal torque wheels 
which receive input from the Earth Sensor Assembly (ESA) for 
EM: and roll and, for yaw, an inertial reference source with 
sun-sensor updates. The ESA is an infrared (IR) sensor that 
views the entire earth and supplies torque input to keep the 
earth centered between four independent detectors. The sun 
sensor uses multiple data inputs from various mechanisms to 
provide the yaw input. 
C. Data Handling Subsystem 

There are four primary components in the data 
handling system onboard NOAA-6; the TIROS Information Processor 
(TIP), the Manipulated Information Rate Processor (MIRP), the 
Digital Tape Recorders (DTR), and the Cross Strap Unit (XSU). 
All the information eventually received on the ground from 
NOAA-6 has to be processed by at least one of these four com- 
ponents. Figure 10 is the data flow diagram for NOAA-6; atten- 
tion is drawn to the path followed by the AVHRR data via the 
MIRP to the switching unit for subsequent transmission at 0.66 
megabits per second to the earth station antenna as real-time 
High Resolution Picture Transmission (HRPT) data. The AVHRR- 
HRPT data stream was the only one used on this project. An 
explanation of the various other sensors shown in the diagram 
can be found in Schwalb (1978). 

The MIRP formats the AVHRR data and adds synchroni- 
Zation, identification, telemetry, and time code information. 
It senses a pulse at the start of each AVHRR scan line and 


initiates a data sampling process that divides the arriving 
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Barth scan data into 2048 computer data words per scan line. 
The pulse that is sensed at the initiation of each scan line 
originates when the AVHRR scan mirror, which rotates at 360 
RPM producing 6 scan lines per second, reaches a precise posi- 
tion in its sweep just prior to scanning across the surface of 
the earth. The data are stored in memory and then subsequently 
read out at a rate suitable for the HRPT on a first-in first- 
out basis. Any one of the 2048 data words or samples, along 
with the number of the scan line on which it is located, de- 
fines a pixel. Throughout this project, the term pixel will 
be defined by the designation (scan line number,sample number) 
or, in short, (NL,NS). A more comprehensive discussion of this 
process can be found in Section III.B below. 
2.  NOAA-6 Onboard Sensors 

There are three primary environmental sensors onboard 
NOAA-6. The TIROS Operational Vertical Sounder (TOVS) con- 
sists of the High Resolution Infrared Radiation Sounder (HIRS/ 
2) whose purpose is to provide data to allow calculation of 
vertical temperature profiles and atmospheric water and ozone 
concentrations, the Stratospheric Sounding Unit (SSU), and the 
Microwave Sounding Unit (MSU). The second sensor, the Space 
Environment Monitor (SEM), consists of a Total Energy Detector 
(TED), the Medium Energy Proton and Electron Detector (MEPED), 
and the High Energy Proton and Alpha Detector (HEPAD). The 
last of the three systems, the AVHRR, was the sensor system 
extensively used on this project and will be described in 


detail below. For an in-depth discussion of the first two 
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sensor systems mentioned above, the reader again is referred 
to Schwalb (1978). 

The AVHRR aboard NOAA-6 is a four-channel scanning 
radiometer that is sensitive to energy in four regions of 
the electromagnetic spectrum. Table 1 below is a summary 


of NOAA-6 channelization. 


Table 1 
NOAA-6 AVHRR Channelization 
(adapted from Schwalb, 1978) 
CHANNEL WAVELENGTH (um) REGION PURPOSE 
1 0.58 - 0.68 visible cloud coverage 


land-water bound. 
snow-ice extent 


2 ئ0‎ 1.1 visible- as above 
near-ir 
p = 3.93 mid-ir sea surface temp. 


cloud mapping 
4 Ое = 11:5 far-ir Sea surface temp. 


cloud mapping 


* 
On NOAA-6, channel 3 is very noisy and usually not used 


An afocal 20.3 cm-aperture telescope, which produces 
EN eld of view of 1.3 + 0.1 milliradians and an instantaneous 
field of view (IFOV) ground resolution of 1.1 kilometers at 
Er (Lauritson, et al., 1979), separates the radiant energy 


ШИЕО СПе four spectral regions with the help of secondary Op” 


tics. The radiant energy in each of these regions is then 
focused on its respective detector. The quantity of energy 
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sensed by the detector then is converted to a count value from 0 
to 255 in the format used for this project. Channel 4 was the 
main channel from which information was gathered and it uses 
a mercury cadmium teluride meto) detector optimized for 
best sensitivity between 10.5 and 11.4 micrometers (Schwalb, 
1978). The spectral response curve for channel 4 is shown in 
Figure 11. The noise equivalent differential temperature 
(NEAT), a measure of the random or coherent two-dimensional 
noise patterns superimposed on the data signal broadcast to 
earth, is less than 0.12 degrees Kelvin at 300 degrees Kelvin. 
Pre-launch AVHRR calibration is covered in a report 
by ITT Aerospace (undated) and post-launch thermal calibration 
of channel 4 is covered extensively in a report by Laurtison, 
шашы (19/9). For every Scan line, the radiometer views deep 
space (0 radiance) and then a blackbody target Gesigned into 
the radiometer housing and kept heated to l5 degrees centigrade. 
To a first order approximation, the radiometer output is linear 
with input energy (Schwalb, 1978) so a two-point linear cali- 
bration, using the above values, is done during every scan 
sequence. Channel 4 with its ао detector, however, has a 
not-quite-linear response due to the physical properties of 
the EE m ron ег а (1979) have generated a table 
of errors for this channel, which represents the difference 
between the actual target temperature and the temperature de- 
rived from the two-point calibration. Table 2 is a summary of 


these data. Note particularly the small errors around 285 


degrees Kelvin, for this is the sea surface temperature range 
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Figure 11.  NOAA-6 AVHRR channel 4 spectral response 
curve (from Kidwell, 1979) 
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determined by the AXBT drops. With this information, a table 
of count-value-to-temperature conversions was generated for 
the time of this project by NOAA-NESS and is included in 


Appendix A. 


Table 2 


NOAA-6 AVHRR Channel 4 Nonlinearity Errors 
Crom Шапты зоп ек а1., 1979) 


TARGET TEMPERATURE ERROR 
Еее К) 
305 U5 
295 0.3 
285 0.0 
275 -0.4 
265 -0.8 


3.  NOAA-6 Orbital Parameters 
Because the development of the satellite data set was 
so closely intertwined with the NOAA-6 orbital parameters, 


their discussion is included in Section III.B.3.c below. 
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met. DATA COLLECTION AND PROCESSING TECHNIQUES 


The use of satellite data on any research topic intro- 
duces extensive data processing problems, especially when one 
considers that a typical NOAA-6 infrared satellite image con- 
tains over nine million pieces of data. This section explains 
the procedures used on this project to collect, process, and 
analyze NOAA-6 satellite imagery with emphasis in the area 
of geographic location accuracy. Also included in this sec- 
tion are the procedures to collect and process the AXBT data 


as well as the collection of the GOSSTCOMP product. 


A. AXBT COLLECTION AND PROCESSING 

As part of the data base for this project, a series of six 
Navy P-3C aircraft flights were staged out of NAS Moffett Field, 
۳۰۰.۶۶113, ×ہ‎ the purpose of dropping a pattern of bathyther- 
mographic sonobuoys (AXBT). The dates of these flights were 
15, 17, and 19 November and 1, 3, and 5 December 1980. These 
flights were scheduled as part of the Naval Postgraduate School's 
research effort on behalf of the joint U.S.-Canadian Storm 
Transfer and Response Experiment (STREX). 

Eact of the nine-hour flights flew northwestward from Cape 
Mendocino, California and proceeded to drop a series of AXBT's 
along a track 1333 kilometers (720 nm) long as shown in Figure 
12. The spacing between the buoys was 55.6 kilometers (30 nm). 
Mme first four flights flew out and back on the center track 


E ihc AXBT's on positions 1 through 24. Flight 5 flew the 
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center track on the outbound leg and flew the southern track 

on the return leg dropping AXBT's on positions 1 through 13 

and 34 through 46. Flight 6 repeated the center track outbound 
E Ev the northern track inbound dropping AXBT's on positions 
rough l3 and 54 through 65. The northern and southern 
return tracks were designed to gather data on the horizontal 
thermal structure and were offset lll kilometers (60 nm) 

either side paralleling the center track. 

The complete navigation suite of the P-3C was used in cal- 
culating the position of each of the deployed AXBT's. At the 
end of each flight, the cumulative error of the inertial navi- 
gation system was checked and recorded. For the flights whose 
data were selected for the project, this error was less than 
4 nautical miles. 

Also important to note is the procedure that the P-3C on- 
board computer uses to calculate the splash points of the de- 
ployed AXBT's. Upon releasing the AXBT from the aircraft, the 
computer ballistics program uses the calculated wind speed 
and wind direction from the navigation system to provide a 
IEectory for the first 2000 feet of fall. After this 2000 
feet of fall, the ballistics program assumes a straight descent 
to the water. This entry point becomes the so-called splash 
point for which geographical coordinates are calculated and 
displayed to the flight crews. Most of the AXBT's were dropped 
from an altitude of 2000 feet except when low clouds or icing 


conditions prevented flying at that altitude. It was felt 





that the location error associated with the few high altitude 
drops was within the 4 nm aircraft navigation error. 
І. The Air-Dropped Expendable Bathythermograph (AXBT) 

The AXBT is an air-dropped expendable bathythermo- 
graph transmitter set deployed by Navy P-3C and S-3A aircraft. 
Its purpose is to provide an accurate profile of the vertical 
thermal structure from the ocean's "surface" to about 350 
meters. Upon water entry, a seawater battery activates, power- 
ing a VHF transmitter, and approximately 30 seconds later a 
temperature probe begins a 5 foot/second descent (Sparton 
Electronics, 1976). The temperature probe and accompanying 
electronics within the sonobuoy package translate the sensed 
water temperature into a frequency broadcast by the radio trans- 


mitter using the formula 
frequency = 800 + 20(temperature deg. F.). 


This low-power broadcast from the sonobuoy is intercepted by 
the aircraft where it is electronically recorded on specially 
processed paper in real-time. 

The accuracy of this process is governed by the accuracy 
of the thermal probe on the sonobuoy and this is claimed to be 
within 1.0 degrees C by the manufacturer (Sparton Electronics, 
1976). Reports in the literature place the repeatable accuracy 
EN Un 0.2 degrees (Barnett et al., 1979). The detailed 
workings of this type of AXBT is described in Sessions and 
Wilson (1976) although the buoy described in their work was 


supplied by a different manufacturer. 


اک 





It is also important to note that the temperature probe 
does not start at the surface of the ocean but begins its 
Есон а depth of about 0.2 meters (Barnett et al., 1979). 
The probe itself also may be subject to very low temperatures 
if the aircraft transports the buoys at a high altitude for 
a long period of time before deployment or if the buoys them- 
selves are dropped from a high altitude. Both Navy aircraft 
described above have systems designed to prevent the freeze-up 
of the sonobuoys, and altitude launch restrictions do apply 
for the deployment of this buoy. 

ША АХБТ Data Processing 

The thermal profiles were recorded using two different 
methods. As described above, the P-3C-produced paper copy of 
the thermal profile was used with a plastic overlay to read 
off the temperature for any depth. These readings were then 
transferred to paper logs by hand. The accuracy of this method 
is within the accuracy limits of the AXBT itself. 

The second method involved the use of an AXBT-digitizer 
provided by the University of Hawaii and which was used also 
in the NORPAX Experiments. This piece of equipment was de- 
signed specifically to be used onboard the P-3C aircraft during 
flights. It connects into the equipment that receives the 
Signal broadcast from the AXBT and digitally records on mag- 
netic tape the signal representing the thermal profile at one 
second intervals. These tapes then are analyzed on a computer 
ious outputs produced. Figure 13 is an example of a 


Group of AXBT profiles produced by this system. Figures 14 and 
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15 show two additional products, a depth-of-isotherm summary 
and a listing of the temperature at 5 meter intervals for 
each AXBT. The accuracy of this system is also within the 
accuracy limits of the AXBT. 

It should be noted also that the AXBT-digitizer had 
provisions for recording the raw received signal from the 
AXBT onto an analog tape recorder. These tapes then were 
used during the analysis phase as a direct input to the 
AXBT-digitizer in order to verify questionable temperature 


profiles. 


SATELLITE DATA SET SELECTION AND PROCESSING 

The decision criteria used to determine which satellite 
passes to examine were reviewed in the following order: 

(1) the satellite pass coverage had to include the ocean 
area where the AXBT's were dropped; 

(2) the time of the satellite pass should be as close as 
possible to the time when the AXBT's were dropped; 

(3) the ocean areas containing the AXBT's should be rela- 
tively cloud-free; and 

(4) there had to be at least one clearly identifiable 
landmark somewhere on the full satellite image. 

Two outside government facilities were used in addition to 
the facilities at the Naval Postgraduate School in order to 
choose satellite passes which met these decision criteria. 

1. NOAA-NESS 

The facilities of the Satellite Field Services Station 


of the National Oceanic and Atmospheric Administration's National 


55 





.15 < ت۶۸7 


date 

























time 


Е. 








Digitizer temperature-at-5-meter-intervals 
summary, an example (from Kilonsky, 


lat. 






























long. 


| 














axbt depth 














no: 















(m) 
| 












1191180 3 192614 423500 1291200 58. 4 жа» 
91481 51413 101483 151403 201491 251483 501481 
301262 551421 601341 691217 701104 751071 801045 
160 933 105 945 110 929 115 917 120 911 125 904 130 902 
139 251 155 043 160 830 165 832 170 827 175 830 189 819 
299 795 205 TBO 210 776 215 773 220 770 225 765 239 760 
250 733 255 “30 260 725 265 717 270 711 275 705 205 698 
309 "074 305 671 310 663 315 662 320 657 325 633 230 657 
5 4 30 260 021 265 622 3T Э 2 
3 1933101 433290 130 800 59 9 
91-24 51445 101442 151446 2921440 251420 201440 
+ 991446 331419 640I311 631136 701091 751059 891049 
100 922 105 959 119 913 115 928 120 912 125 89? 130 900 
150 2525 155 812 160 829 165 822 170 816 175 411 180 893 
229 772 205 со 210 770 215 770 2329 TOU 225 762 239 75% 
A50 723 299 100 209 727 265 122 270 719 275 713 289 706 
209 970 505 671 310 663 315 652 J29 ótt 325 634 339 625 
2509 590 2219 чоо 982 365 504 370 do: 7 
g 19909Ł 442100 131 700 55 6 
01295 51370 10901397 151295 2901295 251395 391397 
501355 95132 601120 651029 TO 991 75 907 09 959 
109 225 105 2C0 110 378 115 BTO 120 892 125 853 130 030 
139 773 155 7065 109 761 105 752 170 744 1735 7238 1609 735 
8E 33 205 ¿35 210 727 215 727 220 727 225 719 239 717 
230 623% 235 632 200 671 265 663 270 663 275 651 289 046 
209 ¿23 303 617 310 606 315 600 220 596 225 592 3139 581 
ВО 271 255 565 240. 561 265 555 370 5 545 
293 49 469 132 700 38 7 
01314 SIDI! 101314 151314 291317 251319 301317 
201217 5515С0 001112 051013 70 988 75 952 H0 933 
156 573 195 851 110 035 115 814 129 794 125 783 139 769 
150 £33 153 793 100 723 165 733 170 793 175 735 169 733 
Bea 205 ¿de 210 727 215 727 220 722 225 722 200 722 
250 701 255 695 260 690 265 001 270 602 275 677 200 6607 
300 $31 303 623 :310 020 315 613 329 608 325 091 120 599 
و رو 3ت2 9ج2‎ 59: 360 330 265 547 JTO 543 > J60 
IS ў 94900 133 800 03 8 
91:70 51500 101306 151390 291396 251303 201202 
3901253 331306 601237 051029 70 960% 75 922 ¿9 906 
109 73% 105 844 119 827 115 7960 129 709 125 742 130 719 
120 053 133 701 160 702 165 702 170 698 175 702 169 705 
209 702 205 (401 2109 695 215 685 220 090 225 687 230 684 
250 653 233 660 7600 052 265 619 270 649 275 641 209 035 
E ^D 565 612 310 СІР 215 (01 320 395 325 599 299 582 
8509:603 350 557 260: 553 265 240 970 542 : TE 
6191145 d 102920 462300 1341200 66 9 
61193 91203 101204 151204 291201 251294 201504 
5612754 33120" 001201 651136 TO 333 75 850 (0 ۳ 
OOO ICS ¿97 110 731 115 769 129 739411 125 742 139 995 
139 922 155 6024 140 701 165 791 170 793 175 703 110 701 
203953 265 600 210 679 215 002 220 602 225 074 229 671 
250 6%9 233 093 260 649 265 6t1 270 639 275 628 289 623 
200 521 305 526 319 390 315 381 320 576 J25 569 339 566 
33. S: 5. есе тот S70 075 : зро 512 





۷1 3 










+6 100 








1334000 64 

91279 91202 1901209 151202 201232 
2391222 251279 001271 051147 170190391 
109 ¿25 105 832 110 019 115 898 129 TOD 
159 722 155 721 160 721 165 524 1750 730 
DKS ECS 210 210 719 215 71% 220 717 
220 637 259 602 260 671 265 666 270 657 
200 029 305 620 310 617 315 609 320 005 
330 374 353 368 360 304 365 557 J79 553 


-— eee н 


- = жу A mo A A 


eras 


56 





282 
913 
2 
тэт 


I< 
637 
601 
545 





2802 
902 
THT 
(27 
711 
652 
293 
042 


301 

9 
130 
16) 
220 
53451) 
$939 


269 


temp. 
(С) 
| 


391482 
051032 


135 
185 
215 
285 
333 


659 
811 
732 
694 
049 


391442 
851034 


135 


163 
235 
285 


335 


094 
TOT 
731 
“00 
014 


351295 


Us 
133 
155 
245 
205 
335 


949 
bas 


* JU 
«903 
өза 
202 


331316 


Go 
135 
165 
233 
235 
335 


920 
TIt 
755 
218 
657 
202 


391396 


55 
195 
155 
245 
2055 


335 


#35 
703 
“06 
619 
631 
974 


391203 


85 
129 
103 
249 
255 


4% 


out 
661 
гол 
667 
617 
299 


351202 


03 
135 
105 
439 
255 


335 


055 
75% 
738 
ot 
өк 
200 


401402 
901013 


140 
190 
240 
290 
340 


801 

797 
740 
688 
644 


401440 
991021 


140 
199 
210 
290 
340 


#55 
7859 
751 
092 
609 


491394 


90 
140 
190 
249 
299 
340 


932 
014 
733 
694 
639 
374 


+01316 


900 
v44 
255 
412 
017 
374 


491305 


90 


379 
0628 
206 
674 
28 
571 


491203 


90 
140 
190 
240 
290 
240 


811 
079 
791 
665 
612 
393 


401201 


99 
140 
190 
249 
290 
J40 


со 
244 
730 
697 
646 
252 


198E) 


77 
451482: 


95 
145 
195 
245 
293 


245 


972 
862. 
794 
740 
682 
644 









77 


451442 


95 
145 
195 
245 
295 
445 


245 


996 
078 
703 
245 
607 
604 









4912314 


95 
145 
195 
215 
295 
345 


089 
735 
733 
710 
639 
971 






451306 


95 
145 
195 
245 
295 
345 


95 
145 
195 
245 
295 
245 


859 
698 
703 
668 
625 
369 





т? 
451204 


005 
604 
698 
663 
606: 
347 








Environmental Satellite Service (NOAA-NESS) in Redwood City, 
California, were used in initially selecting the satellite 
passes. The Redwood City facility is one of three NOAA-NESS 
stations that monitor NOAA-6; the other two are the Command 
and Data Acquisition (CDA) stations in Gilmore Creek, Alaska, 
and Wallops Island, Virginia. Redwood City differs from the 
CDA stations in that Redwood City records the digital High 
Resolution Picture Transmission (HRPT) readout consisting of 
three channels of AVHRR data in the 8-bit precision field- 
station format. These 1600 BPI, 9-track computer-compatible 
magnetic tapes then are archived in Redwood City on a 90-day 
rotating basis. The CDA stations record the various other 
data formats broadcast from NOAA-6 as well as recording the 
HRPT data in 10-bit precision which then are forwarded to the 
NOAA Suitland, Maryland, facility where processing and archiv- 
ing on a more permanent basis occur. The precision loss in 
Eug from the 10-bit HRPT data to the 8-bit HRPT data is 
between 0.4 and 0.5 degrees when making thermal comparisons 
(Kidwell, 1979). The amount of data recorded per satellite 
pass depends on the satellite's elevation in relation to the 
receiving station antenna and can be limited by the l3-minute 
capacity of a standard length magnetic tape. Passing directly 
over Redwood City's antenna, NOAA-6 would be within reception 
Ese for 15.5 minutes, depending on orbital altitude, and 
could provide data from a circular area 6200 kilometers in 
diameter centered on the antenna (Schwalb, 1978). According 


Bewschwalb, the satellite provides useful data only if it is 
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at least five degrees above the horizon. This reduces the 
ӘСЕ time to 13 minutes and the circular area to 5200 
kilometers. 
a. Field-Station Format 

The field-station format differs from the CDA- 
Eon format in that it is a combination ASCII-Binary format 
consisting of a single header record at the beginning of the 
tape followed by up to 15000 data records. Each of the data 


records is a sequential interleavening of the scan lines and 


the recorded channels as shown in Table 3 below: 


Table 3 


Field-Station Format 


RECORD CONTENTS 
1 header 
2 scan line 1--AVHRR channel A 
3 scan line 1--AVHRR channel B 
3 scan line 1--AVHRR channel C 
5 scan line 2--AVHRR channel A 
14998 scan line 5000--AVHRR channel 
14999 scan line 5000--AVHRR channel 
15000 scan line 5000--AVHRR channel 


channel A, B, or C = any sequence of channels l, 2, 3, 4 





The 40-byte header record, all in ASCII, contains 
ашакагоцпа station identification (SFO for Redwood City), the 
channel numbers identifying which three of the four available 
AVHRR channels were recorded, the time (GMT) of the first scan 
line, the duration of the pass, and the orbit number. See 
Figure 16 for an example of the header record. Each of the 
remaining 15000 or so data records have identical 2138-byte 
formats beginning with a 14-byte ASCII "mini-header" consisting 
of an identification sequence, the specific AVHRR channel num- 
ber from which the data in the record originated, the Julian 
date of the scan line, and the time (GMT) of the scan line. 
Following the "mini-header" are 10 bytes of telemetry data, 

6 bytes of back scan data, 10 bytes of space view data, and 
50 bytes of space data, all in binary format. The remaining 
2048 bytes, also in binary, are the video data from which esti- 
mates of the sea-surface temperature are derived. See Figure 
17 for an example of one of these data records. 

b. Ephemeris Data Set 

A set of ephemeris data for NOAA-6 also is main- 
tained at Redwood City. An ephemeris data set consists of 
tracking information so that the field station can capture the 
satellite's data stream as the satellite rises above the hori- 
zon and passes overhead to the opposite horizon. See Figure 
18 and Figure 19 for examples of an ephemeris data set. More 
importantly to this project, the ephemeris also contains the 
subsatellite points for the pass calculated at one minute 


intervals. A subsatellite point is that point on the earth's 
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Figure 16 


Header record--field-station format 


(adapted from Kidwell, 1979) 

WORD BYTE CONTENTS DU TE CONTENTS PYPE 
| iL station ID 2 station ID ASCII 
2 B station ID 4 blank ASCII 
3 5 blank 6 channel A ASCE IL 
+ 7 channel B 8 channel C AS ELL 
5 9 hours 10 hours ACCII 
6 11 minutes 12 minutes ASCII 
7 ES seconds l4 seconds ASCII 
8 155 duration-min 16 duration-min ASCII 
9 17 duration-sec 18 duration-sec  ASCII 

10 ШЭ erbit 20 OLDIE ASCII 

DI 21I s utt: 22 2 ASCII 

l2 23 76 24 blank ASCII 

13-20 25-40 blank 

enenmel А, В, ог С - спаппеі 1, 2, 3, ог 4 


As an example, a pass selected for the project may have a 


header record as follows: 


SFO 1340333481300 7244 


indicating a Redwood City tape (SFO) containing AVHRR channels 


EN апа 4A, Time of the first scan line was 03 hours 33 


minutes and 48 seconds (GMT) while the duration of the pass 


recorded was 13 minutes and 00 seconds. The orbit number was 


7244. 





WORD 
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16-20 
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Figure 17 


Data record--field-station format 
(adapted from Kidwell, 1979) 


БУТЕ CONTENTS BYTE CONTENTS 
1 ID 2 ID 
3 ID 4 П 
5 channel no. 6 day 
7 дау 8 day 
9 hours 10 hours 
Е minutes Ш 2 minutes 
13 seconds 14 seconds 
15-24 telemetry data 
(average) 
2 0ے‎ back scan data 
(average) 
31-40 space view data 
(average) 
41-90 space data 
(raw) 


2۲21۳33٣۱20 


1 


EYEE 

ASCII 
ASCII 
ASCII 
ост 
дос 
ASCII 
ASCII 


Binary 


Binary 


Binary 


Binary 


Binary 





0969:1225 ٤ 
268909969%2%1 
00١86 1 
969و 9ے‎ 56 
909110070 


8٦:96 ١ 


Tete 


و 2599 3339207 


sslNNIn 0+80€£1 * TOL 


ш 3391934 0 7ھ‎ 
= ۸۱۳۰ ۹)۷ 

= 300N “амз257 30 ت۷ ٭‎ 
= ۱۷0٥١۲٢ 141 
* ALIOISJIN3223 


= SIXV VUFVW-IWS3S 


51М54513 1У11840 


! 


°° SWIL H20d3 1V دا‎ ٣۳٢٦٢٢٥٢ NOLANISASY 


ІНІ, Н20643 1Y 


ТЕТЕ.‏ و وج ےج و و جو ےج و و و ےی و یھی 


BVESOSL Lc? C072 
G7S7OR99DLOP ESI 
6799262066” I 
899 64 
182120900070 
COTELESSST°66TL 


NYIu3!d 


096 Я5Л 0351%32 AON Lv 1Y0d38 NOI1VZITVIAINI 


$531*v 0 09 


d 1 1 7 


ATYWGNY NVZW 
359133 30 989 
JOGON 25V düà 357 
0 0 0 0 1 
۸1۱10101003333 


3IXv YOPIW-IW25 


3١ 


3NI1ANOY SSIDIV SIV3N3Hd3 13IVSIZIVAS WVICdA WSLSAS NUILLVOLAYN N-SO811 


(086T 


‘хэхеэла шолз) эташехэ це 


^os елдер ѕтләшәцая 


"gt sıanbrı 


іч 


62 





УЕ ОТ Бе мое 
97 593— S*St— ¿1718 جج٭ہ‎ МЕТ МУ Е $! ? 
м%03- 677 36/1 9%6 носе МА? 2 әр о 
РРО ЭМСЕ сощ 21 М. ТТ 9 9 + 2 6 ۶ 
Lf499—- RI ӨТІЛІ 6222 1 о о 
9$*5G—- 39:22 (УТ (706 MY GIL NOES. 55 st 9 
025-2092 стат аз МНЕ М1%70 56 649 
8” 5° 1a eS cil <a СС MY LITLI N9'657 LS ST 9 
сет 7 ErIE 2°39 9°15 Ме 911 NIS! а GI 7 
0756 1717 CS 2" 9c MI*-lI Уэ”. ц“. І Э 
02687 ҒО СУСЕ Соул Ма а ne Jazz 
сБ «ТЕС УЕЛА ОТ МОСТ NES RASO 
S"t9 L*09 GFZ 9921 Кас Т І. Сб аа; ана? 
5" BO 7250 8090 7557 M6* 801 NyY*95S 15 «І Э 7 
В veto) Me OCT N**^5 Q5 Sl 9 o 
МІТСОІ МРС b7 Sl 9 
МЕ>ООТ №429? 85 61 3 
ML*'9^ MWNM8*^?2 19 $1 چو‎ 
M2* 14 М6%2024 ЭУ ST YD 
MS E. Маа. LY «І 9 
Mot. УИ 3 3 3 ۷۶ 
МОЛ М?Ч Р? 
SIXV  SIXV ЖЕК О 
А X 7 Я A313 YLS tENE Ud. 0S. NI set AVA 
920 39%*97 SNO) 39JN 25% ZLs3/ Sl SI WIL SOON ٣ 15654. *0rOn * A38 


(086T 'xe»xeoag шолу) oetTdwexs ue ‘jas елер этлэшэцая 





"6T әлпБтд 


ike 








surface directly beneath the spacecraft and represents the 
middle of the scan line for the AVHRR. А11 the subsatellite 
points for each scan line taken together represent the ground 
track that the satellite followed in its orbit. It should be 
noted that any alignment errors made when the AVHRR module was 
attached to the spacecraft during construction may result in 
the subsatellite point not being the center of the AVHRR scan 
line. A summary of alignment data may be found in an undated 
report prepared by ITT Aerospace for NASA. For the purposes 
of this project it was decided that any alignment errors were 
so slight as to be negligible and therefore that the sub- 
Satellite point would represent the center of each AVHRR scan 
line. Other information in the ephemeris data set important 
to this project were the orbital elements listed in the preface 
to each ephemeris including orbital period, semi-major axis, 
eccentricity, and inclination. This information was vital to 
me orbital calculations made further on in this project and 
will be explained there. 
C. Initial Satellite Pass Selection 

Each AVHRR scan line is approximately 2840 kilom- 
Es long with 1420 kilometers on each side of the subsatellite 
point. With this information as well as the ephemeris data 
E tor the dates of the P3-C flights and a chart of the 
Northeast Pacific Ocean, it was relatively easy to determine 
specific passes which viewed the ocean areas where the AXBT's 


Were dropped, thus satisfying the first decision criterion 
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mentioned above.  Twenty-two satellite passes were thus 
selected for further screening. 

The selection from these 22 passes of orbits whose 
time matched as closely as possible the time of the AXBT drops 
was done in conjunction with the investigation of cloud cover- 
age over the ocean area of interest. As the project relied 
solely on the use of the infrared channels of the AVHRR and 
because cloud cover effectively prevents AVHRR scan coverage 
of the ocean surface, the absence of cloud cover in the ocean 
area of interest was a major factor in pass selection. Also, 
as will be discussed below, NOAA-6 is a sun-synchronous satellite 
that circles the earth 14.2 times every 24 hours. As a result, 
NOAA-6 views the same earth location at the same local sun 
time each day. This translated to our ocean area of interest 
as between 0330 and 0400 (GMT) for ascending passes and between 
1650 and 1720 (GMT) for descending passes. Since Redwood City 
maintains hourly pictures taken from the visual channels of 
the geostationary GOES-WEST satellite, examination of these 
pictures for cloud coverage in the ocean area of interest re- 
Saeed in the selection of one pass for each of the six flight 
dates that represented the best compromise between matching 
times and cloud coverage. The six passes chosen for further 
examination are listed in Table 4 below. 

2. NASA-Ames Research Center 
ШЫНЫ Il Et Criterion to be satisfied, identification of 
a landmark on each image, was done on the Interactive Digital 


Image Manipulation System (IDIMS) located at the Technology 
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Table 4 
Selected Satellite Passes 
AIRCRAFT PROJ. AXBT БОРЕ ORBIT СЕВЕ ТУРЕ 
FLIGHT NO. TIMES (DTS) NO. TIMES (DTG) PASS 
15 Nov 80 IL 151905-160049 7209 1516539 151 2 D 
Мох 80 2 719322180236 7244 180331-180346 A 
19 Nov 80 3 191843-192232 7266 191712-191716 D 
01 Бес 80 > 011826-020031 7429 010344-010358 A 
03 Dec 80 5 031826-040127 7465 031656=031710 D 
05 Dec 80 6 051816-060013 7486 050355-050408 A 


where D = descending and A = ascending 


Applications Branch of the Airborne Missions and Applications 
Division under the Director of Astronautics, NASA Ames Re- 
Search Center, Moffet Field, California. 
a.  IDIMS 

The IDIMS system is a software package that inter- 
acts with a minicomputer (HP-3000), a display terminal, a 25- 
inch COMTAL display screen, a Dunn Instruments color camera 
recorder, and a high-speed printer, and is used extensively 
to work with satellite data, especially LANDSAT imagery.  Op- 
(5 are available that allow the user to manipulate inter- 
actively satellite imagery so that specific topics of interest 
may be investigated like land use with LANDSAT data, sea surface 
temperature, cloud cover, or ice pack coverage as examples from 
TIROS-NOAA imagery or the many other applications available 


from NIMBUS-7 imagery. The mechanics behind specific options 
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are proprietary data owned by Electromagnetic Systems Labora- 
tory, Inc. (ESL) of Sunnyvale, California, who developed IDIMS 
and to whom the reader is referred for more detailed information 
(ESL, Inc., 1978). A second IDIMS system is located at the 
ВЕЕР Institution of Oceanography in La Jolla, California, 
where initial training was done by the author in the use of the 
IDIMS system. A third IDIMS system, also operated by NASA Ames, 
is located in a mobile van that tours the Western United States; 
its facilities were used upon one of its stops at the Naval 
Postgraduate School. 

b. Landmark Identification 

The basic procedure for landmark identification 
used by this project on the IDIMS system was as follows: 

(1) run a short tape routine on each of the se- 
lected pass's magnetic tapes to identify the number of data 
records, hence the number of scan lines per pass (number of 
data records minus 1 header record, all divided by 3); 

(2) read the 3-channel AVHRR data from the magnetic 
tape into computer memory and initiate IDIMS processing; 

(3) recall that data comprising the infrared 
Channel from memory and display it on the COMTAL display 
screen; 

(4) enhance the displayed infrared imagery for 
temperature using false colors; 

(5) use the Dunn color camera recorder to produce 
ап 8 by 10-inch color Polaroid photograph of the enhanced 


image; 
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(6) use the ZOOM option of IDIMS to enlarge se- 
lected sections of the displayed image in order to locate 
landmark pixels by scan line number and sample number; and 
mally, 

Wise ENES EICERINT operon Of IDIMS to dump to 
the high-speed printer the count values of all pixels within 
a specified area surrounding the landmark. 

An explanation of certain aspects of this procedure 
is explained in the sections below. 

(1) Count Values. In order to display a satellite 
image on the COMTAL display screen, the IDIMS system sequen- 
tially unpacks the video data read into memory from the 
magnetic tape. These data consist of count values between 0 
and 255 which represent the difference in detected energy be- 
tween a look at deep space and a look at a radiating surface 
such as the earth or clouds (Schwalb, 1978). These count 
values are used by the IDIMS system to produce an image with 
a grey-scale intensity range from 0 to 255 in order to match 
the same range of count values. Options available on the IDIMS 
System allow various color assignments based on these count 
values including an automatic full-spectrum false color assign- 
ment where red is "hot" and blue is "cold" or vice versa. Also 
available are options allowing a single color to vary in its 
Saturation over the full count range of 256 values or the 
assignment of a specific color to an individual count as 
Sepue=count 125" or to a range of counts as "blue=counts 125 


BuEough 130". Assigning colors in this manner tends to produce 
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a confusing image if many hues are used indiscriminantly. 
For purposes of this project, the automatic full-spectrum 
false color assignment of red "cold" and blue "hot" was used 
so that the oceans were blue and the cloud tops, being much 
colder, were red in all the photographs. 
€. Pixel Identification 

The identification of a landmark pixel and the 
subsequent assignment of a scan line number and sample number 
are made easy by the IDIMS system but the principles behind 
their assignment had to be understood so that other landmarks 
and buoy positions could be located as needed later on in the 
project. 

As seen in Table 4 above, three of the six selected 
NOAA-6 passes were ascending passes and three were descending 
passes. An ascending pass is one where the satellite in its 
orbit crosses the earth's equator heading northwards while a 
descending pass is one where the satellite crosses the equator 
heading southwards. As the satellite is moving, the AVHRR 
scan mirror sweeps from right to left perpendicularly across 
the satellite's velocity vector six times per second with each 
sweep defining a single scan line. Although the scanning 
mirror rotates in a complete circle, only the data located 
55.4 degrees either side of nadir is retained. Nadir is a term 
Similar in meaning to the subsatellite point in that nadir 
represents that point on the sweep of the scanning mirror 
when the mirror is pointed at the spot on the earth's surface 


directly underneath the spacecraft. The radiometer data 
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stream from this 110.8-degree arc is electronically divided 

by the MIRP into 2048 equal samples each representing 0.054 
degrees of the total arc (110.8 divided by 2048). Therefore, 
when the satellite transmits its data stream to the earth 
receiving station, the first 2048 bytes of video data taken 
together is intrinsically labeled scan line number 1 while the 
first byte is intrinsically labeled sample number 1 ог, as 

used throughout this project, (1,1). The second byte of the 
first scan line is labeled (1,2). Sample number 1505 from 

scan line number 3520 would be labeled (3520,1505). Sample 
number 1 through 1024 are located to the right of the sub- 
satellite point when looking in the direction of the satellite's 
velocity vector while samples 1025 through 2048 are to the left. 
Because of this right-to-left pixel numbering system, when an 
image is displayed on the COMTAL unit by the IDIMS system, an 
ascending pass image looks reversed and upside down while a 
descending pass image looks normal where normal is defined as 
having Alaska to the north or top of the image, Hawaii to the 
west or left of the image, and California to the east or right 
of the image. An ascending pass image, when displayed on a 
display screen or when stored into computer memory for further 
processing, has Alaska on the bottom, Hawaii on the right, and 
California on the left of the image. This occurs because IDIMS 
always displays pixel (1,1) in the upper left corner of the 
COMTAL unit. Although landmark identification can be easily 
done on either type of pass, the geometry involved further on 


in this project rests heavily on a clear understanding of which 
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type of pass you are analyzing and on which side of the sub- 
satellite point is the landmark or buoy. Figure 20 is an 
example of an ascending pass while Figure 21 is an example of 
a descending pass. In each figure, the satellite would travel 
directly up or down the center of the image respectively. 

The IDIMS system automatically keeps track of this 
numbering system and conveniently displays the scan line num- 
ber, sample number, and count value for the pixel you have 
identified on the COMTAL unit using a movable cursor. By use 
of the ZOOM feature on IDIMS and with reference to a chart of 
the local area, it is usually easy to identify landmarks. In 
most of the passes used on this project, the San Francisco Bay 
Area was identified readily and its (scan line number, sample 
number) determined. For purposes described later on, up to 
20 landmarks per satellite pass along the west coast of the 
United States from Glacier Bay, Alaska, south to Mexico and 
from San Francisco east to Pyramid Lake, Nevada, were identi- 
fied. Only one of these landmarks is needed to navigate the 
image as will be explained later on. 

The PICPRINT feature of IDIMS allowed the dumping 
to a high-speed printer of the count values of the pixels sur- 
rounding the landmark. This was done as a method of verifying 
the accurate position of the feature chosen for landmark iden- 
tification. Figure 22 is an example of a PICPRINT output where 
by using a variation of the game of connecting the dots one 
can connect count values in order to recognize features. This 


EEnod of verification will not work, or is made more difficult, 
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NOAA-6 ascending pass, 


Figure 20. 
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UCP RINE output of count values, 
cloud-free image 


Francisco Bay entrance; 


Gure 22. 
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if the temperature of the land is the same as the water, if 
i lying clouds (fog) have similar temperatures as either 
Е water or land, or if clouds obscure sections of the land. 
milly, the visual channels of the spacecraft are used but 
since many of the passes occurred at night, only the infrared 
channels were usable. Figure 23 is an example of a Case 
where clouds interfered with the landmark identification 
process. 

Because an average pass contained 4680 scan lines 
Beer with 2048 samples or about 9.6 million pixels per pass, 
a convenient method of matching pixel numbers (hence count 
values) to AXBT positions was necessary; thus a system was 
needed to "navigate" the image. 

Satellite Image Navigation 
E-wWentioned Sn the introduction, there are many sources 

Of error when one wants to compare satellite-derived sea sur- 
face temperatures and AXBT-derived sea surface temperatures. 
It was decided at the beginning of this project that an at- 
tempt would be made to reduce as much as possible one of 
these, that being the earth location errors associated with 
transferring AXBT drop positions to a satellite image, so that 
temperature comparisons could be made. Several methods were 
tried and three of those methods that produced the smallest 
errors are described below. For the purposes of developing 
procedures for locating AXBT's on the satellite image only, 
an assumption was made that the geographical location of the 


AXBT's was accurate; hence any aircraft navigation errors, 
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Figure 23.  PICPRINT output of count values, San 
Francisco Bay entrance; cloud-covered 
image. (Note: This image is from an 
ascending pass hence it appears upside 
down.) 
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ballistic errors on the falling AXBT's once launched from the 
TE, drift errors on the floating AXBT, or human errors 
in transcribing positional data from aircraft displays to logs 
were ignored. These sources of error will be discussed later. 
a. Zoom Transfer Scope 

A Bausch and Lomb Zoom Transfer Scope was used 
meeerally in an attempt to transfer the AXBT positions to 
the satellite image. A zoom transfer scope allows one opti- 
Sally to overlay a chart, on which the AXBT positions have 
been plotted, onto a satellite image where enough distinguish- 
ing features (landmarks) are evident so that by optically 
۳٣۰166 7 condensing, or rotating the chart, landmarks on 
both chart and image coincide. Once the landmarks coincide, 
the operator manually marks with a pencil the AXBT positions 
onto the satellite image. While the system works fine with 
small area images consisting mostly of land, it could not be 
used satisfactorily on this project for a number of reasons. 
First, each of the selected NOAA-6 passes covered an area ex- 
tending from Northern Mexico to Alaska and from mid-Pacific 
ENDS western United States. Reducing a chart of this area 
ЕС а size suitable for use on a zoom transfer scope (about 10 
by 10 inches) necessarily requires reduction in the accuracy 
m lotting geographical coordinates. Second, on each of the 
NOAA-6 passes, approximately 80 to 90 percent of the coverage 
area was open ocean with any visible landmass only on the edge 
of the image.  Landmasses on the edges of these images are 


much more distorted than landmasses near the subsatellite 
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Ems due to a combination of satellite scan geometry, earth 
curvature, and the transfer of these images to a flat medium 
like a photograph or chart. Third, the AXBT's were dropped 
along a line over 1300 kilometers long stretching northwestard 
from Cape Mendocino, California. There are no landmarks in 
the Northeast Pacific Ocean between Cape Mendocino and the 
Aleutians; therefore location accuracy decreased the farther 
away from the coast the AXBT's were dropped. Fourth, marking 
a chart manually with a pencil necessarily involves inaccura- 
cies especially when one is trying to locate geographically 
an item as small as an AXBT. Last, and the hardest to over- 
come, is that once the buoy position is marked on the satellite 
image, some method must be found to determine the pixel number 
and hence the count values of the AXBT's position.  Remember- 
ing that there are 9.6 million pixels per image, determining 
the exact pixel to choose for a count value would involve some 
guesswork and possibly even large-scale pixel averaging. 
Satellite images unfortunately do not come marked with latitude 
and longitudes, nor do charts contain scan line numbers and 
sample numbers. 

0۶ ۰ 9 4 

A second method of trving to locate an AXBT on 

the satellite image involved the use of an ESL, Inc-developed 
IDIMS function called TRNSFORM.  TRNSFORM is used mainly in 
registering LANDSAT imagery and involves the calculation of a 
transformation matrix between matching sets of control points 


Using a least-squares fit method (ESL, Inc., 1978). A first, 
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second, or third order transformation is possible.  TRNSFORM 
was not designed to navigate NOAA-6 imagery mainly because 
TRNSFORM requires 15 to 20 landmarks spread over the entire 
image in order to obtain a small pixel error. Therefore, use 
of this function also failed to provide the accuracy desired 
for this project for one of the same reasons that the zoom 
transfer scope failed, in that landmasses were present only 
on the edges of the images. 

The method finally used, and from which a location 
accuracy of less than 2 pixels resulted, was the development 
of a computer program that determined a satellite's orbit 
referenced to a single landmark in the image and from this, 
when given an AXBT latitude and longitude, could determine the 
scan line number and sample number of the AXBT. The development 
of the program required a basic understanding of the orbital 
dynamics of NOAA-6 as well as a working knowledge of spheri- 
cal geometry. 

C.  NOAA-6 Orbital Dynamics 

Кот Ога! information in ehis section, the work 
by Stewart (1979) and Schwalb (1978) was used extensively. 

The NOAA-6 satellite is a sun-synchronous satellite 
which means that its orbital plane rotates at the same rate 
as the rotation of the earth about the sun. As a result, the 
satellite views a point on the Earth's surface at the same 
local sun time each day. Table 4 above listed those times 
that NOAA-6 viewed the AXBT drop area. According to Schwalb 


(1978), the orbital plane precession rate is approximately 
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КОМСОО 000000199 radians per second or 0.986 degrees per 
day eastwards. This rate is achieved by placing the satellite 
Ші ап огріс with a suitable inclination. Іп the case of the 
NOAA-6 satellite, the inclination was determined prior to 
launch to be 98.739 s 0.15 degrees, where inclination (i) is 
defined as the angle the satellite's orbital plane makes with 
the earth's equatorial plane measured counterclockwise from 
east. А retrograde inclination (i ) is the supplement of 


the inclination. See Figure 24. 


Figure 24 


NOAA-6 orbital plane inclination 
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The period of the satellite, obtainable from the 
ephemeris data set (as is the inclination), is the amount of 


Е it takes the satellite to make one orbit of the Earth. 


30 


The predetermined launch value for the NOAA-6 period was 101.58 
minutes; therefore, NOAA-6 orbits the Earth 14.18 times per 24 
hours. For each orbit the earth rotates 25.40 degrees east- 


ward. See Table 5 for a summary of NOAA-6 orbital parameters. 


Table 5 


NOAA-6 Orbital Parameters 


orbital plane precession rate 0.986 deg/day east 
Шісісппасіоп (1) 98.739 = 0.15 deg 
Beerograce inclination (1 ) O ol SOS deg 
period 101.58 minutes 
orbits per day 14.18 

earth rotation per orbit 25.40 degrees east 
orbital altitude о ап 


The predetermined launch altitude of the satellite 
ELE. - 18.5 kilometers. The orbit of NOAA-6, to a first 
approximation, is an ellipse. From the ephemeris data set 
the semi-major axis of the ellipse and its eccentricity can 
be found, thus making the satellite's altitude on any pass 
simple to calculate as shown in Figure 25 below. The mean 
Satellite altitude (H) is that distance measured from the 
center of the Earth and can be derived using Equation (1). 
Using the above orbital information, the data in 
Table 6 below was derived for the six selected NOAA-6 passes 


ШЕЯ in this project. 
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Satellite Altitude Determination 








a semi-major axis (from ephemeris) in nm 
e eccentricity (from ephemeris) 
Р И 2152 
Ы semi-minor axis з а[(1-е )] / 
Н mean satellite altitude = (2b+a)/3 (Eq. 1) 
Table 6 
Satellite Data Set Orbital Parameters 

ORBIT 7209 7244 7266 7429 7465 7486 
PARAMETER 
inclination ( degrees) 98.69708 98.69708 98.69708 98.69123 98.69123 98.69123 
retrograde inclination 81.30292 81.30292 81.30292 81.30877 81.30877 81.30877 
period (minutes) 101.13285 101.13285 101.13285 101.13084 101.13084 101.13084 
semi-major axis (km) 7185.4875 7185.4875 7185.4875 7199.1856 7199.1856 7199.1856 
eccentricity 0.001187 0.001187 0.001187 0.000603 0.000603 0.000603 
mean satellite altitude 7185.4841 7185.4841 7185.4841 7199.1847 7199.1847 7199.1847 
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The reason that the orbital parameters are not 
constant for each pass is that the satellite is subject to 
Eu forces that tend to cause its orbit to vary. The largest 
of these forces is the fact that the earth is not a perfect 
sphere but an oblate spheroid.  King-Hele (1958) and Brouwer 
(1959) developed mathematical solutions to describe this per- 
turbation whose primary effects on the orbit include changing 
the orbital plane precession rate and changing the period. A 
secondary influencing factor is the effect of atmospheric 
drag on the satellite which acts to change the eccentricity 
and is a function of the satellite's altitude. A third influ- 
ence is the effect of solar wind and radiation. It should 
be noted that during 1980, the International Solar Maximum 
Year, solar flare and sunspot activity reached some of the 
highest levels recorded (Ponte, 1981). Lesser influences in- 
clude the gravitational effects of the sun and the moon on 
the satellite. 

d. Computer Navigation Program 

The main computer navigation program was developed 
under the following premise: given the orbital parameters of 
NOAA-6, the latitude and longitude of an AXBT, and a satellite 
image upon which one landmark has been identified as to (scan 
line number, sample number) and latitude and longitude, calcu- 
late the (scan line number, sample number) of the AXBT so that 
the count value, hence the sea surface temperature, of that 


pixel can be identified readily either on IDIMS or any other 
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computer system.  Procedural methods were outlined by Mueller 
№0981). 

(1) Preliminary Programs. Three preliminary com- 
puter programs were designed to be run on the IBM 3033. The 
first program, SCANLINE, was a simple block counter that counted 
the number of data records on each magnetic tape. The number 
of data records minus the header record divided by 3 gives the 
total number of scan lines per pass. The program listing can 
be found in Appendix B. The second program, TAPEDUMP, was a 
routine designed to dump from the magnetic tape any number of 
bytes per data record and to translate their ASCII-Binary for- 
mats into decimal notation. This program was used to verify 
that there were indeed six scan lines per second and its 
listing can be found in Appendix C. The third program, AREAMAP, 
was designed to function in a manner similar to the IDIMS' 
PICPRINT function in that it would extract from the magnetic 
tape the count values of a selected grouping of pixels around 
ВЕ landmark or AXBT pixel. This program was used to verify 
landmark locations and to determine the surface thermal struc- 
ture around the position of the AXBT. Its listing can be 
found in Appendix D. 

(2) Common Case Geometry. In the development of 
the main computer program, it was necessary to consider four 
Gases in the process of predicting an AXBT pixel. These four 
cases are: 

(a) an ascending pass where the landmark 


has a sample number greater than 1024 (Case 1); 
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(Б) an ascending pass where the landmark 
has a sample number less than or equal to 1024 (Case 2); 

(с) a descending pass where the landmark 
has a sample number greater than 1024 (Case 3); and 

(d) a descending pass where the landmark 
has a sample number less than or equal to 1024 (Case 41). 

Common to each of these four cases was the 
assumption of a spherical earth with a radius equal to the 
earth's radius at the landmark latitude. This local radius 


can be calculated using Equation (2) 


2 2 
(3443.925) ++ ٦ ; 
where: 
R = local earth radius in nm 
Ба = landmark latitude in degrees. 


Also common to all four cases were the calcu- 
lations to determine the great circle distance between the 
subsatellite point and the pixel containing the landmark. 
These calculations refer to Figure 26 below. 

These calculations are made possible under 
the assumptions that the subsatellite point is directly be- 
EN the satellite on the earth's surface, that the scanning 
mirror of the AVHRR forms a scan line perpendicular to the 
Beeellite's velocity vector, and that the earth is a perfect 


sphere. 
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Figure 26 


Determination of great circle distances 








H 


subsatellite point 


landmark 
pixel 


where: Н = mean satellite altitude from Eq. 1 


R = local earth radius from Eq. 2 


8 = scan angle 

= great circle distance 
= geocentric angle 

= zenith angle 


Determination of the scan angle (8.) in degrees 
assumes an equal division of the arc viewed by the radiometer 


(110.8 degrees) into 2048 samples, tnus 


: (landmark sample-1024)(110.8/2048) if sample > 1024 
S = 


(1025-landmark sample) (110.8/2048) if sample < 1024 


The zenith angle (9) in degrees now can be found: 


ren Fe 
p R 
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The geocentric angle m in degrees is found from Equation 


)3 ےت و - .8 = 8 


If desired, the geocentric angle in degrees can be expressed 


as the great circle distance in nautical miles from Equation 


(4). 
$ = +8 ٦ (Eq. 4) 


(3 Units and Notation. Because the IBM 3033 
and the Fortran computer language were used heavily during 
this project, all angles were converted to or used in radians. 
Table 7 lists the common conversion factors used.  Notations 
on all figures included in this project were designed to have 
the same definition whenever possible so comparisons could be 


made between the four cases. 


Table 7 


Program Conversions 


45 degrees = I radians = (8) [(tan (1.0 radians)] 
T 
any angle in radians - (same angie in degrees) (,/4) 


any geocentric angle in 
radians expressed as a 
great circle distance 
in nautical miles 


= > (geocentric angle) (60) 


NE uticalemilo ^ 1.835 kilometers 
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(4) Nodal and Subsatellite Point Calculations. 

The main program, LOCATE, is divided into two sections. The 
first section calculates the orbital characteristics referenced 
to the previously-identified landmark. The desired output is 
the time and longitude of the ascending or descending node. 
Their calculation is dependent on the four cases enumerated 
above and described in detail below. Once the time and longi- 
tude are known, the second part of the program can proceed to 
calculate the pixel number for an AXBT. To begin, the calcu- 
lation of the ascending or descending node and time follows 
for each of the four cases. 

Case l. Ascending pass with landmark sample number greater 

than 1024. 

The derivation of orbital characteristics 
referenced to a single landmark in Case 1 made use of Figure 
27 below. 

E ime CGAL 0 пе Satellite is directly 
over the subsatellite point that has an unknown latitude (Lo) 
and longitude A). The only known quantities are that the 
scan line that includes the subsatellite point also includes 
the landmark with known latitude (LI; longitude (43); and 
from IDIMS, a known scan line number (NL); and sample number 
INS). From the ephemeris data set for this pass, the inclina- 
tion (and hence retrograde inclination (i_)) and the period 
are known. From Equation (3) or (4), the great circle dis- 
tance "s is known. The goal of this orbital set of calcula- 


tions is to find the latitude and longitude of the subsatellite 
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Figure 27 
@asewl—-oroltal characteristics 
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point, the longitude of the ascending node и) and the 
time of the ascending node. 
By using similar triangles and the Law of 


Sines, the angle (s) can be determined as follows: 


sin 


: € 
from triangle I AZ وت‎ 
sin L sin p 
O O 
| sin (ІІ. -с ДЕ 
arom triangle II зіп (1. - £) = ; 
sin 0 sin 6 
O 
-1 Sin l= 
hence, = = tan 
sin o 
E + COS 1 
sin Lo = 
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From triangle I 


sin L 
к ГО 


From triangle II, the distance over which the satellite 
travelled between the ascending node and the subsatellite 


point is 


COS 
?o 


% سکس‎ [cos $ 
g 


Finally the latitude of the subsatellite point (2) can be 
determined using Equation (5) and triangle III 


ib = sin”? [sin(i )sin($,)] . (E 3) 


Por the moment, the rotation of the Earth is ignored, so the 
change in longitude at the equator between the landmark longi- 
tude A) and the fixed ascending node longitude (А5 --where 
the "f" indicates a fixed or non-rotating earth derived term) 
can be found by solving triangle I as follows: 
COS $ 
a Е z! O 
و‎ К сова. 
O 

The change in longitude at the equator between the subsatellite 
point and the fixed ascending node longitude can also be found 


from triangle III 





Now, the longitude of the subsatellite point I9 can be 


calculated easily using Equation (6) 


EsNSccount for the earth's rotation during the time the 
satellite traveled from the ascending node to the subsatellite 
point, the longitude change due to rotation (AA) is calculated 
and subtracted from the fixed ascending node longitude as 
follows: 


ф 
At (seconds) = x= (period in seconds) , EO) 


= 


where At is the time the satellite took to travel between the 


ascending node and the subsatellite point. Continuing, 
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/ = — 2 A 
AA 51. (1.002738) (16) (Eq. 8) 


where 1.002738 is the sidereal day correction factor. Finally, 





The time of the ascending node can be found by subtracting the 
At in seconds from the landmark time in seconds. Landmark 
time can be read from the magnetic tape using the tape dump 
program mentioned above to dump the data record containing 
the scan line of the landmark. 
Case 2. Ascending pass with landmark sample number less 
than or equal to 1024. 

The orbital calculations for this case are 
very similar to those of Case l and have the same goals. Figure 
28 shows the geometry applicable in this case. As in Case 1, 
the known values are the great circle distance m from Equa- 
mens 3 or 4 and the landmark's latitude (Lo), longitude D 


and the scan line number and sample number (NL,NS). 


Figure 28 
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Using spherical triangles I and II, the angle 


(е) can be determined as follows: 


Dom triangle I sin (teat) = zur ; 
SIN L sin ¢ 
O O 
ОТОП triangle IL ERE ur E Е ; 
sin ?g sin De 


therefore, equating the two equations and solving for e yields 


| san i 
= tan”? : = 
Y sin L 
sin ےپ‎ = COS l. 
g 
Semeinuing with triangle II 
| _ біті біп 9 | 
о sin Е : 


Ehe distance the satellite travels from the ascending 
node longitude to the subsatellite point of the landmark's 
scan line IF also can be found by 


tanto 
= cos [ 


Е 


(Sin с) (сап фо 


From spherical triangle III, the latitude of the subsatellite 
point can be found using Equation 10 
-1 


L. „л ц 7 (sin ШЕЛ 3 (EG 10) 


2 








The change in longitude between the longitude of the ascend- 
ing node and the longitude of the subsatellite point now can 
m@eetound from triangle III, ignoring the earth's rotation. 


tan L 
S 


] 


о ance 


The angle (a) from triangle III is 


4 f 
ES Sin ۵۸ ے‎ 


| n 
a = Sin [ 


sin Ф. 


] 


which can be used in triangle IV to find the change in longi- 


tude between the subsatellite point and the landmark (AA): 


(cos a) (sin q ) 


cos L 
O 


The longitude of the subsatellite point (А2) now is found 


anng Equation ll 


AN (Dem, 


A = A Е (БО. 22) 
Bene Equations (7) and (8) from Case 1 to determine the 
@eqrees of longitude through which the earth turns during the 


meme the satellite travels from the ascending node to the 
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subsatellite point (AA), the ascending node longitude e) 
Ере found using Equation (9). The time cf the ascending 
node is found exactly as in Сазе І. 
Case 3. Descending pass with the landmark sample number 
greater than 1024. 

The case for the descending pass is slightly 
different from the ascending pass cases in that the satellite 
Mmemtravelling from north to south. Тһе goals and the known 
factors are identical with the ascending pass cases. Figure 
29 below pertains to the spherical geometry applicable in this 


case. Note that the ascending node becomes the descending 


node (dn) in descending pass calculations. 


Figure 29 
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The calculations for this case are exactly 
the same as those for Case l with some exceptions as noted 
below. The fixed earth change in longitude between the land- 
mark longitude and the descending node longitude Gn ayy can 
be determined from 


1,505 Фо 


= cos [————] . (e en 


AX an cos L 
O 


The change in longitude between the subsatellite point longi- 
tude and the descending pass longitude (AA) can be found from 
۲۰۰٦1916 1 


COS 6$ 
-1 
AA = cos [ E 


cos L 
S 
The longitude of the subsatellite point 8 can be determined 


from Equation 14 


n E E 
О E шм), (Eq. 14) 


and the fixed earth descending node longitude now also can be 


determined 


Using Equations (7) and (8), the degrees of longitude through 
which the earth turns while the satellite travels between the 
subsatellite point and the descending node can be determined 


(ДА), and from this the rotating earth descending node 
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longitude can be found using Equation (15) 
A = A FAA аш) 


The time of the descending node now can be determined by 
adding the time calculated in Equation (7) to the landmark 
time. 

Case 4. Descending pass with landmark sample number less 

than or equal to 1024. 
In this final case, the goals and the known 

quantities are the same as in the other cases described above. 
Figure 30 is used to describe the geometry associated with 


this case. 


Figure 30 
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The calculations for this case are exactly 
the same as those for Case 2 with some exceptions as noted 
below. The distance the satellite travels from the subsatellite 


point to the descending node, $ becomes 


i 


COS 0 
E -1 O 
RES los 7 


] 


The latitude of the subsatellite point, Lo, now can be found 


using Equation (16) as follows: 


г. = sin ![(sin i.) (sin 821. (Ба. 16) 


The fixed-earth change in longitude between the landmark 


longitude and the descending node longitude, ie now can 


be found by 


ff a No 
Man = “os [ 


І, 


cos L 
O 


and the change in longitude between the subsatellite point 
longitude and the descending pass longitude, Aà, can be found 


Tom triangle III 


COS ф 
-1 
AX = cos | 5 


cos L 
$ 


The longitude for the subsatellite point, Е is now found 


using Equation 17 





and the fixed-earth descending node longitude becomes 


f 5 E 
A n E ان‎ 1 Е 
ИП using Equations (7) and (8), earth rotation is con- 


sidered now, with the rotating earth descending node longitude 
described by Equation (18) 

аа = Хар + AA à (КО. 719) 
The time of the descending pass is determined by adding the 
time calculated in Equation 7 to the landmark time. 

Ориол ра ел шаёдеЕ са Езоп. With the calcu- 
lation of the time and longitude of the ascending or descend- 
ing node, the second part of the program can proceed to calculate 
MNE, NS) of the AXBT. The previous set of calculations 
referenced the orbital characteristics to the landmark pixel 
(remember pixel = (NL,NS)). The geographical relationship 
between the latitudes and longitudes of the landmark and the 
AXBT are known so the purpose of the remaining part of the 
program is to transform this geographical relationship into 
Satellite image coordinates of (NL,NS). Besides the quanti- 
ties calculated above, the only other known quantity is that 
the AXBT has a unique (NL,NS). 

As a first guess, any arbitrary scan line can 
be chosen to be the "true" scan line containing the AXBT. One 


of the 2048 samples along the "true" scan line could be the 


29 





sample number of the AXBT. The aim of the calculations below 
is to prove or disprove, geometrically, that the arbitrary 
scan line is the true AXBT scan line and, once the true scan 
line is selected correctly, to calculate the correct sample 
number. 

The procedure to determine the time that the 
satellite recorded the scan line containing the landmark was 
described above. Since the time of the ascending or descend- 
ing node was calculated in the earlier part of the main pro- 
gram, subtracting the two times describes the satellite flight 
time between the particular node and the subsatellite point of 
the landmark scan line. The assumption was made earlier that 
the satellite's orbit can be considered circular with a mean 
altitude (H) and that the period of the satellite was the 
amount of time it takes the satellite to complete one orbit; 
then the flight time between the node and the landmark scan 
line can be described as a great circle distance in radians by 


4 _ Рае time) | (Eq. 19) 


S period 


Another assumption was made earlier, which was proved by using 
one of the preliminary computer programs described above, that 
the satellite records six scan lines per second. Combining 
these factors, it is simple to determine the difference in 

scan lines between the landmark scan line and the first-guess, 
arbitrarily-chosen scan line; divide by 6 to get the time 
difference between the two scan lines; either add to (ascending 


EES) or subtract (descending passes) this time difference 
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from the flight time; and use Equation (19) to determine the 
distance between the particular node and the subsatellite point 
of the arbitrarily-chosen scan line. The latitude and longi- 
tude of this subsatellite point now can be determined using 
spherical triangles. 
Assuming that the arbitrary scan line is very 
close to the true AXBT scan line, the next step is to find 
the correct sample number along this line. Beginning with 
sample number 1, and then taking every 89th sample (1,90,179, 
...,1959,2048) the distance between the subsatellite point 
and the center of the pixel can be found using Equations (1) 
through (4) for each of the 24 sample numbers. With this 
distance it is possible to calculate the latitude and longitude 
of the center of each of these 24 samples by using one of the 
four cases described below. 
Case l. Ascending pass with sample number greater than 1024. 
This case uses the geometry of Figure 3l. The 


angle (е) can be found easily by 


_] 5 p 
€ = COS ae | ; (БО. 20) 
S 
Brom which 
і, = cos "[(cos L4) (sin К (EG. 21) 


The distance, in radians, of з. сап be found by using 





۶۲۰16 1ٰ 


Case l--pixel determination 
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^p As 
1 sin L. 
ra einen T Б ; (542 20) 


hence the latitude of the possible AXBT position n can 


be calculated 


L, o= sin "[(sin (p,-+,)) (sin і 0] 


The angles (y) and (АА) now can be determined by 





and 


-і (біп y) (sin 9,) 
AA, = | | == 


Ges № 
S 
thus the longitude of the possible AXBT position L is 


determined simply by 


Case 2. Ascending pass with sample number less than or 
equal to 1024. 


Case 2 geometry uses Figure 32. 


Figure 32 
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The angles (e) and (1) and the distance ШЕ) 
are found using Equations (20), (21), and (22) respectively. 


The latitude of the possible AXBT position becomes 


5 mul. EM ES 
Lo = SIN [sin (¢ ۾‎ + VL sin(i.)] 


The angle (AA P) now is found by 


(SIN E) (Sin фа! 


J1 , 


COS L 
Р р 


and therefore the longitude of the possible position is 


Case 3. Descending pass with sample number greater than 
1024. 
Басс Шс зц лс Figure 33 for its 
geometry, uses exactly the same equations as in Case 1 with. 
the exception that the longitude of the possible AXBT position 


ms round by 
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Figure 33 
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Case 4. Descending pass with sample number less than or 
equal to 1024. 
Case 4, although Figure 34 represents its 
geometry, uses exactly the same equations as in Case 2 with 
the exception that the longitude of the possible AXBT posi- 


tion is found by 


mm GGG со calculating latitude and 
longitude for each of these 24 samples along the arbitrary 


Ene; it also is necessary to calculate the great circle 
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Figure 34 


Case 4--pixel determination 
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distance (d) between the AXBT geographical position and those 


of the 24 samples using Equations 23 


= | 1 1 سج‎ » 
d = cos ¡(sin L, sin L_J+ (cos O MSN gl, (КООШ 2) 


where: 


Г, ene ue AXBT latitude 


23) 


ti 
к 


true AXBT longitude . ( 


By taking the sample number that has tne 


emaliest of these great circle distances, creating a bracket 
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89 samples wide on either side of this center sample, and pro- 
ceeding as before through the appropriate case number for each 
of the 180 samples in this bracket, the sample that has the 
shortest great circle distance between itself and the AXBT 
position can be selected. The reason for the wide bracket 
is to allow for earth rotation and curvature whose effects are 
especially noticeable on the edges of the satellite image. If, 
as before, we assume that we had a scan line very close to 

the true AXBT scan line, a 10 by 10 pixel "Square" box is 
created around the sample with the smallest great circle dis- 
this box is not perfectly square due to 


Bee, In reality, 


the curvature of the earth and the motion of the satellite 
during the scan sequence. Those boxes near the subsatellite 
point would be more perfectly "square" than those boxes on the 


edges of the image. 


The scan lines on 
box as well as the samples on each 
subjected to the same calculations 
the latitude and longitude of each 


corners. The geometry of this box 


The calculations to 


the top and bottom of the 
Greene four corners are 
described above to find 

of the samples on the four 
is shown in Figure 35. 

ов ME 


find the (NL,NS) 


AXBT begins with finding the change in the latitudes and 


longitudes with respect to the changes in the sample and scan 


iue numbers. 
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Figure 35 
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м = longitude of sample on i-th corner 
Lori = latitude and longitude of the sample 
2 having the smallest great circle dist. 
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Continuing, two intermediate equations are required 


кі » 3۸ эл | 

А = An Ay) + 21587 (sample А) + Ey (scan line x); 
Бы 2L (sample A can ine x) 

Pag 9 (NS) 3 (NL) ” 


Finally, determination of the (NL,NS) of the AXBT can be 


made 
A _ В 
۰ 94/3 (NS) JE WS) 
31/3(NL) _ 3L/3 (NL) : 
3۸/3 (NS) 3L/3 (NS) 
A _ r3۸/3 (NL) 
NS = 3۸/3 (NS) 7775085) on 


Rarely will the arbitrary scan line chosen 
as a first guess be close to the true scan line of the AXBT. 
In this case, after the great circle distances have been calcu- 
lated between the initial 24 samples and the AXBT position, and 
the "square" box has been set up, a check is made to see how 
"small" is the smallest great circle distance. As described 
earlier, the satellite follows a ground track as it travels 
poleward that cuts across meridians of latitude at an angle sec 
ШУ 165 inclination. This means that scan lines are not 
oriented east-west along degrees of longitude but are oriented 
northwest-southeast (descending pass) or northeast-southwest 
(ascending pass) crossing many degrees of longitude and lati- 


tude. The check involves comparing the latitude of the AXBT 
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with the latitudes Lo and L, as shown on Figure 35 above. If 


d 
the arbitrary scan line is very close to the true AXBT scan 
line, for ascending passes, latitude Lo will be south of the 


EET latitude while latitude L, will be north of the AXBT 


d 
latitude. The reverse is true for descending passes. If 

the arbitrary scan line is far away from the true AXBT scan 
line, both La and Ба will be either north or south of the AXBT 
latitude. In this case, the smallest great circle distance 

is converted to an integral number of scan lines which are 
added to or subtracted from the first-guess arbitrary scan 
line number depending on whether L_ and L, were Boch 8 ۶2 
north of the AXBT's latitude respectively. The jump to a new 
scan line initiates the entire procedure again beginning with 
the steps necessary to calculate Equation (19). This jump 
process terminates when the number of scan lines to be jumped 
is 5 or less at which time the boxing procedure begins with 
the eventual determination of the AXBT's (NL,NS). 

The main program, whose listing may be found 
in Appendix E, was initially set up to be run interactively on 
a display terminal. The program was used to locate all the 
AXBT's that were dropped from the P-3C. Verification of the 
accuracy of the program was done by selecting 15 to 20 land- 
marks per image and asking the program to predict the (NL,NS) 
even though they were known already from the IDIMS system. 


Results of this verification will be discussed below in Section 


Ш.А. 
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eee GOSS TCOMP 

The GOSSTCOMP sea surface temperature charts were obtained 
for the period of this project from NOAA-NESS. These charts 
are produced on a weekly basis by NOAA-NESS using procedures 
 Епес by Brower et al., (1976) and since updated to take 


advantage of the AVHRR on NOAA-6. 
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ID RESULTES 


A. NAVIGATION ACCURACY 

A major effort was made on this project in an attempt to 
reduce the effects of geometric distortion associated with 
locating landmarks or open-ocean positions on satellite imagery. 
Previously published works with earth location errors in ex- 
cess of 10 kilometers at nadir were suitable for regional 
location and analysis of mesoscale features; however, it was 
believed that accurate comparisons of thermal data were signi- 
ficant only if the products being compared were co-located in 
the same geographical position. The location of thermal 
features is especially important in naval tactical applications. 
If a submarine were taking advantage of the unique acoustical 
properties of an eddy or ocean front, then acoustical prosecu- 
tion by the opposing forces would be more successful if the 
sensors employed by this group were located so as to take 
advantage of the thermal feature also. If pre-mission informa- 
tion mislocated the edge of the front or eddv due to the 
geometric distortion inherent in satellite imagery, then the 
results could be disastrous to one of the parties. Ultimately, 
any error in sensor placement could prove disastrous whether 
caused by satellite imagery or not; part of the purpose of 
this project was to make the location error as small as possible. 


Using the main computer program, the (NL,NS) of each AXBT 


was predicted. The error in this prediction was determined 
Beebe within 2 pixels of the true AXBT (NL,NS). The procedure 
2 





to verify this accuracy began with the identification of the 
(NL,NS) of up to 20 landmarks on each satellite image. Each 
of the landmarks then was treated like an AXBT and its geo- 
graphical coordinates were input to the computer to see what 
the program would predict for each landmark's (NL,NS). These 
predictions then were compared to the IDIMS-determined 
(NL,NS), and the separations in pixels were determined. The 


results are summarized in Table 8 below. 


Table 8 


Statistical Summary of Navigation Accuracy 


SCAN LINE ERROR 


mean 1.32 scan lines 
standard deviation 0.93 scan lines 
992 confidence level 0.88-1.90 scan lines 


SAMPLE ERROR 


mean 1.35 samples 
standard deviation l.17 samples 
99$ confidence level 0.31-1.90 samples 


The conclusion drawn from this statistical summary was 
Maat the predicted AXBT (NL,NS) is within 2 pixels of the true 
AXBT (NL,NS). Because of earth curvature, pixels close to 
nadir are not as wide as those out on the edges of the image. 
Sample number 1 and 2048, on the right and left edge of the 
image respectively, are 4.3 kilometers wide whereas sample 


number 1024 and 1025, located to the right and left of nadir 
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respectively, are only 0.77 kilometers wide. As a result, 

the 2-pixel error can be as small as 1.9 kilometers, if the 
EBEdijcted (NL,NS) is at nadir; or as large as 10.7 kilometers, 
if the predicted (NL,NS) is at the edge of the image. Table 

9 lists the navigation error associated with selected sample 
numbers. Notice that the error is not linear with distance 
from nadir but is less than 5 kilometers over 80$ of the 

image and less than 3 kilometers over 50$ of the image. Only 
on the outer 10% of the image does the error balloon from 4.8 


to 10.7 kilometers. 


Table 9 


Navigation Errors Associated with a 2-Pixel Error 


SAMPLE NUMBER ERROR (km) 
1 007 
200 48 
400 0 
600 5 
800 20 
1024 1,459 
1200 2520 
1400 ЖБ 
1600 30 
1800 4.3 
2048 3007 
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There are some other sources of navigation error that 
should be kept in mind when using data developed by this method. 
Alignment of the AVHRR module during its assembly prior to 
launch could be the source of constant ofíset error. This 
type of error was described previously in this paper; analy- 
sis ОЁ the navigation results did not show any consistent off- 
Bias that could be attributed to module alignment errors. 

Hough the process of verifying the 2-pixel accuracy, 
many landmarks were identified on the IDIMS system as explained 
above. The data from NOAA-6 infrared channel number 4 were 
meee tor landmark identification and their use could introduce 
errors in assignment of the (NL,NS). These errors arise from 
trying to identify: landmarks whose surface temperature may 
not be very different from the surrounding surfaces. This 
effect would become even more pronounced if ground fog were 
present. When selecting these landmarks, the best contrast 
was effected by land-water boundaries, examples of which were 
Point Lobos west of Monterey, the San Francisco Bay entrance, 
Alcatraz Island, Point Reyes, the Columbia River mouth, Lake 
Tahoe, and Glacier Bay among others. Many possible landmarks 
were not considered if there were insufficient contrast to 
ify the feature. A good example of this was the Seattle- 
Tacoma-Olympia area where the numerous bays and tributaries 
had surface temperatures close to land temperatures, thus 
Making it very difficult to distinguish a specific pixel as 


being some peninsula or promontory. 
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A third source of error involved the number of signifi- 
۴۰۰۰٠٠۰۶۰ gures used in the mathematics of this project. Single- 
bit precision was used during computer processing. Although 
this may have had an effect after numerous computations (the 
average program run to calculate 24 AXBT positions per image 
executed 475,000 statements), it was felt that the number of 
significant figures was more critical. An example of this was 
the determination of decimal geographical coordinates. The 
navigation system on the P-3C supplied the computer-calculated 
coordinates of the AXBT's to seconds of latitude or longitude. 
One second of latitude error is equal to 0.1 kilometers, 
which in itself is not so large; however, most landmarks were 
meeneitied using charts with scales of 1:2,000,000. After 
determination of the coordinates, a decimal conversion to three 
decimal places was completed. If errors in this procedure 
were compounded by weak land-water contrast on the infrared 
image used for landmark identification, it could contribute 
Significantly to the 2-pixel error. 

A fourth source of error has to do with the resolution of 
the AVHRR itself. As discussed earlier, the 1.1 kilometer 
resolution would necessarily make it difficult to identify 
something like the Transamerica Building in downtown San 
Francisco. An example of where this could be a contributing 
factor to the 2-pixel error would be in using the most western 
point of Point Reyes. If the scan sequence is such that the 
radiometer does not resolve this point, then the first pixel 


it does identify as being land would be to the east of the 
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BEL The user of the satellite image would have a very 
difficult time in trying to determine whether or not this has 
occurred. As a result, the user would assign the geographi- 
cal coordinates of the most western point, introducing a l- 
pixel error immediately before any computer processing begins. 
It is felt that if landmark's could be identified using sharper 
land-water contrast or using a visible channel vice the infra- 
БЕП Channel if it is available, and if geographical coordinates 
could be assigned with greater accuracy, the majority of the 
2-pixel bias would be eliminated. 

Another error to be considered is that the satellite may 
Abe perfectly stable in its orbit. It is likely that small 
ВЕ Of pitch, roll, or yaw occur from time to time although 
the ADACS system was designed to keep these attitudes to a 
minimum. 

The last error to be discussed is the use of several assump- 
tions made during this project. The Earth was assumed to be 
spherical and although the radius was calculated to be that 
Eus at the latitude of the landmark, a small error will be 
EeEoduced in calculations involving the earth radius term at 
the latitude of the buoys. Similar small errors arise with 
the assumptions made concerning the satellite orbit during 
Feen line calculations, and with the calculation of the mean 
altitude of the satellite above the earth's surface. 

in conclusion, it is believed that the 2-pixel error found 
on this project could be reduced further to sub-pixel accura- 


cies if some of the errors described above were eliminated 





or refined. Since the development of LOCATE, a program using 
most of the same techniques as LOCATE has been developed to 
predict the geographical coordinates of open ocean images 


with similar accuracies (Mueller, 1981). 


B. THERMAL COMPARISONS 
MM Horizontal Distribution 

As usually can be expected when satellite-derived sea 
surface temperatures that are uncorrected for atmospheric 
attenuation are compared with AXBT values of sea surface 
temperature, the satellite temperatures were colder than the 
AXBT data by a mean difference of 2.9 degrees C. Table 10 
lists the mean temperature difference values and the corres- 
ponding standard deviations for this and following comparisons. 
ums 36, 37, 38, and 39 show sea surface temperature com- 
parisons of this and other methods to be described below along 
Ehe buoy line. The majority of the 2.9-degree error can be 
attributed to the effects of the intervening atmosphere be- 
tween the ocean's surface and the satellite radiometer. Cloud 
contamination of the satellite values was not considered a 
major factor due to the screening process that went into se- 
lecting the data. Out of the six satellite passes selected 
for study at the beginning of this project, only three met the 
full requirements that were required for processing. Two of 
the passes were not considered due to scattered clouds over 
enough buoy positions to make any comparisons useless and one 


pass was not considered because it contained no clearly 





Table 10 


Temperature Comparison Statistics 


COMPARISONS 
MEAN (C) STANDARD 
satellite vs. AXBT DEVIATION 
17 November -3.0 05 
01 December -2.6 Шен 
05 December -3.0 0-6 
overall -2.9 On 
satellite vs. GOSSTCOMP 
17 November -2.0 0.9 
01 December -3.2 0.6 
05 December -4.3 0.9 
overall -3.2 0 
GOSSTCOMP vs. AXBT 
17 November -0.7 О 
01 December 0.5 0. 6 
05 December JE 9:26 
overall 0c 3 4770 
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Figure 36. Sea surface temperature comparisons, 
17 November 1980, center track 
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Figure 37. Sea surface temperature comparisons, 
] December 1980, center track 
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identifiable landmarks. Passes selected for complete pro- 
cessing were 17 November, l and 5 December 1980. From personal 
observations onboard the P-3C during AXBT deployment, it was 
noted also that there was no ground fog to interfere with 
satellite measurements of the AXBT positions. The time 
difference of one, three, and three hours between the last 
buoy drop and the satellite flyover for 17 November, 1 апа 

5 December, respectively, probably is not a factor as the 
lowest and highest mean error values were found on 1 and 5 
December with 17 November having an intermediate value. If 
there were a correlation, one would expect 17 November to have 
the smallest error but this was not the case. The transient 
warming of the surface waters during the afternoon, the so- 
called afternoon effect, did not occur due to the weather 
conditions during the three-week project period; hence, this 
process also was ruled out as a source of the error.  Con- 
stant wind speeds in excess of 20 knots from the south on 17 
November, in excess of 25 knots from the northeast on l Decem- 
ber, and in excess of 15 knots from the westnorthwest on 5 
December (National Weather Service, 1980a) along the buoy 

line kept the surface waters under constant wind-mixing condi- 
pons. in addition to the winds, ship observations of the 

sea state at the northern end of the buoy pattern £ound four 
to ten foot swell and two to six foot waves. These turbulent 
mixing conditions are diametrically opposed to the formation 


of the afternoon effect (James, 1966). 
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A method of "field-calibrating" the satellite data to 
eliminate the effects of the atmosphere was suggested by 
Tabata and Gower (1980). Using a simple linear regression 
technique, they plotted numerous ship-obtained surface tempera- 
tures versus satellite count values and found that over a 
limited area and a limited time period between satellite and 
ship observations (1.5 days), the error between satellite and 
ШЕ values could be reduced to 0.5 degrees C. This tech- 
nique was tried using the data from this project. An absolute 
mean difference of 0.3 degrees (s = 0.2 degrees) was found 
between satellite and AXBT values. The time period between 
Satellite and AXBT observations was three to eight hours. 

DMUIShteesews usually will not have the luxury of ex- 
Bending 24 AXBT's on a tactical mission however, so the linear 
regression technique was tried using only two buoys. The 
rationale behind using two buoys was that this is the number 
ОЕ AXBT's usually carried on both S-3A and P-3 aircraft. Addi- 
tionally, a scenario could exist whereby a satellite photo 
mea prior to the flight could pinpoint two sections of 
the tactical operating area where thermal differences exist 
and those two locations could be designated for AXBT deployment. 
The point is to try and get a spread in temperature between 
the two AXBT's. Using the two-buoy method, the data from AXBT 
positions 1 and 13 on 1 December were used for the linear re- 
gression. Predictions of temperatures from count values found 
mean difference error of 0.3 degrees (s = 0.3 degrees), the 


same value found using the 24-buoy method. When this same 
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regression formula with constants calculated from 1 December 
data was used to predict 5 December (4 days later) and 17 
November (15 days earlier) temperatures, mean errors of 0.5 
degrees (s = 0.3 degrees) and 0.4 degrees (s = 0.3 degrees), 
respectively, resulted. Although examination of all possible 
cases would be necessary before conclusive results could be 
stated, these preliminary estimates indicate that it should 
be possible to use two AXBT's and an infrared satellite image 
to predict sea surface temperatures within 0.5 degrees for at 
least two days after the original satellite pass, a prediction 
tool particularly helpful if clouds obscure the sea surface 
during those two days or if AXBT assets are in short supply. 
In addition, these procedures can be used in near real-time 
processing of current satellite images and do not rely on any 
atmospheric model processing. In any case, these thermal 
predictions that are very accurate in location and fairly 
accurate in temperature would be tactically significant, in 
place of mean values or best-guess values, when doing sound 
velocity calculations near meander, eddy, or frontal regions. 

Relative temperature gradient analysis displayed the 
expected correlation between satellite and AXBT values. Both 
the AXBT and satellite gradients were 0.6 degrees per 60 nm 
on 17 November and 1 December while on 5 December the AXBT 
gradient was 0.56 degrees per 60 nm and the satellite gradient 
was 0.52 degrees per 60 nm. 

An interesting thermal feature whose horizontal sur- 


face manifestation was detected by both the satellite and 


126 





the AXBT was a meander in the final stages of closing off 
from its parent body of water to form an eddy. See the darker 
region transected by the buoy line through buoy positions 23, 
3, 22, and 4 in Figure 40. This warm-core meander had an 
approximate 100 nm diameter with the exception of an open arm 
extending southward into its parental water mass. The diameter 
was verified by the buoys dropped on the northern and southern 
tracks. These buoys were 60 nm away from the center rack and 
the warm meander did not show up on any of the thermal traces. 
The satellite indication of this diameter resulted in a slightly 
larger radius, a fact attributed to the thermal resolution limi- 
tations of the satellite data in determining weak temperature 
boundaries. In the satellite images, this meander is surrounded 
on the west, north, and east by the Subarctic Current-California 
Current confluence. The center of the meander had a surface 
temperature of 15.8, 14.9, and 14.4 degrees C on 17 November, 
І and 5 December, respectively. A chart of the monthly mean 
surface temperature for November 1980 (Renner, 1981) clearly 
shows the intrusion of a large tongue of warm water from the 
area between San Francisco and Hawaii northward along the 
west coast of the United States. 

The decrease in the surface temperature of the meander 
over the project time period was reflected by the decrease 
in the surface temperatures all along the buoy line. Both 
the satellite and the AXBT's recorded mean changes of 0.9 
degrees between 17 November and 1 December and 0.7 degrees 


between 1 December and 5 December. This drop in temperature 
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is to be expected considering the weather conditions as men- 
tioned above. Оп 17 November, a low pressure system was 
firmly entrenched over the Aleutians while a high pressure 
system was anchored off of Southern California. This pressure 
pattern is typical of the Northeast Pacific in early winter. 

A cold front extending southward from the Aleutian low moved 
across the buoy line during the evening of 17 November and 
crossed the U.S. coastline during the morning of 18 November. 
Winds before passage were southerly at 25 knots while after 
passage the wind shifted northerly at 30 knots; hence the 
condition for considerable wind-mixing existed. A series of 
cold fronts on 20-22 November and 25 November also passed 
through the project area continuing to lower the sea surface 
temperature and drive the mixed layer depth deeper. Оп 1 
December, low pressure cells were established west of the 
coast of Washington and about 500 nm west of the central 
California coast. The Washington low strengthened and cen- 
tered near buoy positions 7 on 2 December. This low was 
accompanied by winds in excess of 35 knots on 3 December over 
the entire buoy area while a slow-moving cold front hugged 

the coastline. On 4 December, a high pressure area, previously 
metab lished in the Gulf of Alaska, moved into the project 

М О trom the north pushing the cold front well inland although 
mae) low remained off the Washington coast. The high moved 
southerly on 5 December, influencing the weather over the 


entire buoy area (National Weather Service, 1980b). See 
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Figures 41, 42, and 43 for the surface weather depiction 
charts for 17 November, 1 and 5 December, respectively. 
Because the linear regression model mentioned earlier 
was used to remove the effects of the atmosphere with fairly 
good results, a comparison was made between the satellite 
data and the GOSSTCOMP product. Sea surface temperature pro- 
ducts from GOSSTCOMP have been subjected to an atmospheric 
correction model and are issued on a weekly basis. On all 
three days, the satellite data from this project were colder 
than the GOSSTCOMP values. The mean difference for 17 November, 
1 апа 5 December were 2.0 degrees (s = 0.9 degrees), 3.2 
degrees (s = 0.6 degrees), and 4.3 degrees (s - 0.9 degrees) 
respectively with the overall mean of 3.2 degrees (s = 1.1 
degrees). This overall mean agrees fairly well with the 3.5 
to 3.9 degree bias enumerated by Klein (1979). The reason for 
the 3.2 degree bias can be attributed directly to the effects 
of the atmosphere, exactly the same situation as seen in the 
AXBT versus satellite comparisons mentioned previously. An 
interesting point to be made is that the 3.2 degree bias of 
GOSSTCOMP versus satellite data is higher than the 2.9 degree 
bias of AXBT versus satellite data. This led to a comparison 
between GOSSTCOMP and AXBT data with the result that GOSSTCOMP 
values were 0.3 degrees (s = 1.04 degrees) warmer overall than 
the AXBT values. Because the project area was never totally 
cloud-free during the period of observations, 1t 1s felt that 
the overcorrection for atmospheric effects described by Klein 


(1979) is still a factor in the warmer GOSSTCOMP values. PE 
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should be mentioned that the GOSSTCOMP product did not 
indicate the warm meander that the AXBT and satellite data 
Meedecad, probably due to the large grid structure used by 
GOSSTCOMP. 

All three methods of sea surface temperature determin- 
ation, AXBT, satellite, and GOSSTCOMP, were compared to the 
20-year mean surface temperature values of Robinson (1976). 
AXBT values versus climatology resulted in the November AXBT's 
being 0.3 degrees colder than the mean while the December 
values were about the same as the mean. The reason for the 
cooler surface waters in November is probably a result of 
the high incidence of weather frontal passage with accompany- 
ing high winds through the project area. As was expected, 
climatology did not show the warm meander. 

Comparison of satellite versus monthly mean data 
found the satellite data averaging 3.0 to 2.7 degrees colder 
than the mean for November and December. Comparison of GOSSTCOMP 
versus the mean resulted in GOSSTCOMP being 0.4 degrees warmer 
than the mean for November and 2.0 degrees colder than the 
mean for December. 

Ze Vertical Distribution 

There is no known way at present to sense remotely 
the vertical thermal structure in the ocean; however, if one 
combines knowledge of the horizontal gradients with clima- 

БКО Осу, а fairly accurate synopsis of the upper ocean thermal 


structure is possible. A more accurate picture can be 
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formulated if the satellite data are augmented with well- 
placed AXBT drops. 

From climatology, the expected mean surface tempera- 
ture and the mean layer depth for November were 12.9 degrees 
and 50 meters (s = 5 meters) respectively over the project 
area. The AXBT mean surface temperature and mean layer depth 
for 17 November were 12.6 degrees and 58 meters (s = 6 meters). 
For December, climatology means were 10.8 degrees and 67 meters 
(s = 6 meters) and the AXBT means were 10.8 degrees and 71 
meters (s = 7 meters). Figures 44, 45, 46, and 47 show the 
vertical structure along the buoy line on 17 November, 1 and 
5 December (center track and north track) respectively. From 
the numerous oceanographic studies in the area (Tully, 1961; 

NE 1961; etc.), it is known that during this period of 
the year, the layer depth is deepening towards the maximum 
limit of the top of the permanent halocline at 100 meters. 
The 100-meter depth is not reached usually until February. 
The deepening of the layer is directly attributable to the 
turbulent mixing conditions caused by the sustained high wind 
speeds and by the convective mixing caused by the surface 
cooling during the calmer periods. The weather pattern for 
late-November and early-December was discussed previously. 

A general rule of thumb is that warm surface waters generally 
exhibit shallow layer depths while colder surface waters ex- 
hibit deeper laver depths. This pattern held true throughout 
the project area. Although definitive sea surface temperature- 


layer depth relationships were not within the scove of this 
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ЮЕСЧЦЧЕСС, it is conceivable that flight crews, after observing 
a horizontal and vertical thermal analysis done on AXBT data 
and satellite images from a previous day(s) could view cur- 
rent satellite infrared images and through comparisons of 
surface temperatures (atmospherically corrected or not) could 
estimate the layer depth. 

Further deployment of layer depth pređiction tech- 
niques may be found in combining the works of James (1966), 
McAlister and McLeish (1970), and Mollo-Christensen and 
Mascarenhas (1979). James described a wav of using heat bud- 
get calculations and wind-mixing parameters to predict the 
ocean thermal structure.  McAlister and McLeish described an 
airborne system that was capable of measurement of the total 
heat flow from the sea. Mollo-Christiansen and Mascarenhas 
used LANDSAT data to calculate heat storage in the upper mixed 
layer of the ocean. The use of satellites to estimate wind 
speed and direction has been demonstrated. If the heat budget 
could be estimated using the principles described in the papers 
above, and the winds determined, then the procedures described 
by James may be applicable. 

An example of how present satellites are inadequate to 
mse the vertical structure is the fact that only the AXBT 
traces located a region of a subsurface temperature maximum 
on the northern end of the buoy line. This region would help 
to define the lower extent of a subsurface sound channel. 


There was no surface manifestation of this feature, like the 
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warm meander, and hence the satellite completely missed 1t. 
This temperature maximum was relatively narrow (50 meters) 
and had an axis between 150 and 175 meters. On 17 November, 
the area affected by the temperature maximum was relatively 
large with the axis at the shallower depth. By 5 December, 
the area affected was reduced by half with the axis depth 
deepening in conjunction with the layer depth. The area of 
the subsurface temperature maximum can be seen in Figures 44, 


fey 46, and 47. 
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V. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


It was determined that geographical locations over open 
ocean areas could be located on satellite imagery to within 
2 pixels or within 3 kilometers over 50% of the image (within 
5 kilometers over 80$, worst case within 10.7 kilometers). 
Location of satellite features on geographical charts using 
a Variation of the main program has similar accuracies. 

Satellite observations can be a very effective tool if 
used in conjunction with available groundtruth data. The 
development of the NOAA-6 AVHRR allows a more accurate de- 
termination of sea surface temperatures than from previous 
satellites; however, the satellite mean values differ from 
groundtruth values by 2.9 degrees C. This bias is attribu- 
table directly to the effects of the atmosphere. By using a 
method of linear regression, it is possible to "field-calibrate" 
the satellite data so that the mean error between satellite 
and AXBT data is reduced to 0.3 degrees. This small bias 
held true whether 24 or 2 buoys were used to define the cali- 
bration constants for the linear regression equation for a 
single satellite pass. When these same constants were used 
to predict the sea surface temperature on satellite passes 4 
days later and 15 days previous, the mean error increased 
slightly to between 0.4 and 0.5 degrees. This three-week pre- 
diction period saw the passage of numerous cold fronts and 


experienced long periods of turbulent mixing conditions so the 


142 





relatively small bias may make this procedure a possible 
ШКО сп tool in calibrating satellite images without the 
need to use atmospheric attenuation models. 

Relative temperature gradients were constant, as ex- 
pected, between satellite and AXBT data values. The intro- 
duction into the project area of a warm-core meander in the 
final stages of becoming an eddy was apparent readily from 
both the AXBT data and the satellite imagery, although it was 
ignored by GOSSTCOMP. The detection of subsurface thermal 
structures by satellites is successful only if there is a 
surface manifestation of the feature such as the meander. A 
subsurface temperature maximum that is associated with a sub- 
surface sound channel was not detected by the satellite, as no 
surface manifestations were present. 

When GOSSTCOMP values were compared to AXBT values, an 
overall 0.3 degree bias was noted, although individual daily 
values ranged from 0.7 degrees colder to 1.2 degrees warmer 
than the corresponding AXBT data. 

Over the project period, the sea surface temperature 
decreased on the average 1.8 degrees while the layer depth 
deepened on the average 13 meters. The cooling and deepening 
process was directly related to the number of frontal passages 
and the high wind speeds during the project period. It may 
be possible to use systems currently under development to sense 
remotely the existing thermal structure to 100 meters; or to 


use a combination of satellite instruments to determine the 
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heat budget and wind speeds and from this information to 
calculate the thermal structure. 

Determining the exact location of thermal features on 
satellite imagery is only the beginning in tactically employ- 
ing satellite data for naval missions. Further development 
of the relationship between the sea surface temperature and 
the vertical thermal structure is warranted. An excellent 
starting point would be along lines similar to those developed 
under the ASWEPS program. If a possible relationship could 
be derived, the use of environmental satellites by naval 
tacticians could advance far beyond its present stage. 

Using methods developed by this thesis, it is possible to 
develop hypothetical scenarios that need to be tested opera- 
tionally.  Supposing that the 17 November satellite imagery 
and AXBT data, and the l December satellite imagery only were 
available, an antisubmarine aircraft flight crew assigned a 
surveillance mission on 1 December could have predicted within 
one half degree the sea surface temperature over a wide ocean 
area. Knowing both the vertical thermal structure from 17 
November and the effects of cooler surfaces and sustained 
fem winds on this structure, the flight crew could make a 
fairly accurate prediction of the 1 December oceanographic 
and acoustic conditions. Development of techniques that relate 
the surface temperature to the vertical structure would make 
this process that much easier. 

2۰۱۱٠۹ co predicting the vertical structure, some 


aspects of the horizontal structure are equally important, 
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especially to the fast-paced antisubmarine warfare efforts 

of the carrier-based S-3A aircraft.  Submarines may be able 

to use sharply-delineated fronts and eddies to keep themselves 
acoustically hidden from aircraft carriers while remaining 
within weapons firing range. The persistence of these thermal 
features over a period of time is seen easily in satellite 
images and their exact geographical location can be determined 
using the methods derived in this thesis. A current satellite 
photo is far superior to other products now available in con- 
veying this type of information. As an example, GOSSTCOMP 
missed totally the small (100 nm diameter) warm-core meander 
Sound п this project.  GOSSTCOMP is also not a real-time 
Product. 

In conclusion, the following recommendations for further 
study are suggested: 

(1) the development, where applicable, of an empirical 
relationship between the sea surface temperature and the 
vertical thermal structure; 

(2) the development of thermal structure prediction tech- 
niques using both satellite data and the empirical relation- 
ship developed above; 

(3) a study of the persistence of horizontal thermal 
features using satellite imagery; 

(4) a study of how accurately satellite-observed surface 
thermal features reflect the subsurface structure; 

(5) the development of a streamlined LOCATE program 
suitable for ship-board use so that surface thermal features 


from satellite imagery can be located more accurately; 
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(6) the further development and testing of the "field- 
calibration" technique of dealing with atmospheric attenua- 
Баса: апа, 

(7) the investigation of including a current satellite 
image, with thermal features geographically located, in 


mission planning information. 
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APPENDIX A 
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